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L  INTRODUCTION 


Empirical  thermospheric  models  are  employed  for  a  number  of  practical  uses 
including  satellite  ephemeris  prediction.  For  users  who  are  unfamiliar  with  the  field  of 
thermospheric  modeling  it  can  be  difficult  to  choose  among  the  many  models  which  h.«ve 
been  developed.  This  report  is  intended  to  help  such  a  user  select  the  model  best  suited 
for  his  requirements  by  presenting  a  summary  of  the  characteristics  of  a  number  of  the 
most  important  empirical  models  of  the  thermosphere.  With  a  view  toward  use  in 
satellite  orbit  studies,  the  primary  emphasis  is  placed  on  total  mass  density  rather  than 
temperature  or  number  densities  of  individual  species.  An  earlier  report  by  Bruce  (1965) 
summarized  empirical  models  in  use  as  of  1965. 

The  next  section  of  this  report  will  describe  the  most  important  empirical  thermo¬ 
spheric  models.  A  synopsis  of  early  models  is  followed  by  a  more  complete  description  of 
those  in  current  use.  These  descriptions  will  be  followed  by  a  section  which  provides  a 
comparison  and  discussion  of  the  relative  advantages  and  limitations  of  the  models  in 
current  use.  Finally,  a  comprehensive  set  of  figures  which  illustrate  the  model 
predictions  of  total  density  as  a  function  of  solar  and  geophysical  parameters  at  several 
altitudes  of  interest  will  be  given.  These  figures  will  allow  a  more  detailed  comparison 
of  specific  quantitative  aspects  of  the  various  models. 

The  effects  of  aerodynamic  drag  on  satellite  orbits  are  most  important  in  the 
altitude  region  from  120  to  400  km  above  the  surface  of  the  earth.  Satellites  which  have 
perigee  altitudes  below  120  km  suffer  such  great  deceleration  that  their  lifetimes  before 
reentry  are  so  short  as  to  be  impractical.  For  satellites  with  perigee  altitudes  above 
400  km,  the  orbital  perturbations  caused  by  drag  are  relatively  small,  and  the  resulting 
gradual  change  in  orbital  elements  can  be  accounted  for  easily. 
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The  atmosphere  above  120  km  is  characterized  by  a  spectacular  increase  in 
temperature  with  increasing  altitude,  as  seen  in  Fig.  1,  and  it  is  this  feature  from  which 
the  region  derives  its  name,  the  "thermosphere."  At  high  altitude  the  temperature 
finally  levels  off,  approaching  a  limiting  value  known  as  the  "exospheric  temperature."  A 
number  of  models  have  been  developed  in  an  effort  to  represent  the  behavior  of  the 
thermosphere.  These  models  generally  fall  into  one  of  two  classes,  theoretical  or 
empirical.  The  theoretical  models  attempt  to  achieve  a  description  of  thermospheric 
behavior  in  terms  of  the  basic  physical  processes.  Although  it  is  through  such  models 
that  an  understanding  of  thermospheric  dynamics  must  be  obtained,  for  many  uses  these 
theoretical  models  are  not  suitable  due  to  their  complexity  and  the  large  amounts  of 
computer  memory  and  time  they  require. 

Empirical  models  of  the  thermosphere  make  little  or  no  attempt  to  describe 
correctly  the  underlying  physics,  but  rather  try  to  provide  density  values  in  a  convenient 
manner.  It  is  these  empirical  models  which  are  most  suitable  for  applications  such  as  use 
in  operational  systems  concerned  with  the  effects  of  aerodynamic  drag  on  satellite 
orbits. 

Prior  to  1957  there  had  been  a  limited  number  of  rocket  probes  of  the  thermo¬ 
sphere,  but  there  were  insufficient  data  to  determine  the  systematic  density  variations. 
An  Air  Research  and  Defense  Command  (ARDC)  model  was  developed  in  1956  which 
utilized  the  data  from  these  probes.  However,  only  the  average  density  and  temperature 
were  presented.  Analysis  of  the  decay  of  satellite  orbits  due  to  aerodynamic  drag 
allowed  study  of  the  variations  of  atmospheric  density  near  perigee  and  it  was  soon 
confirmed  that  there  were  substantial  variations  of  density  at  satellite  altitudes.  The 
major  variations  in  density  are  due  to  changes  in  the  amount  of  heating  of  the  upper 
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atmosphere  by  solar  extreme  ultraviolet  (EUV)  radiation.  Although  near  120  km  the 
variations  in  density  are  small,  at  400  km  the  density  may  change  by  over  two  orders  of 
magnitude  from  the  minimum  to  the  maximum  of  the  eleven  year  solar  cycle.  TTie 
simple  ARDC  1956  model  was  therefore  of  obviously  limited  value.  A  number  of  models 
based  on  satellite  data  were  later  developed  to  describe  the  observed  variations  in  an 
orderly  manner. 

CIRA  1965  was  the  first  widely  used  model  which  included  all  known  major 
thermospheric  variations  modeled  continuously.  The  diurnal  variation  and  variations  due 
to  changes  in  the  level  of  solar  activity  (as  characterized  by  the  10.7-cm  solar  flux)  were 
included.  Also  incorporated  was  the  heating  of  the  high  latitude  thermosphere  during 
geomagnetic  storms  by  energetic  particles  which  causes  increases  in  the  density,  some¬ 
times  by  a  large  factor  in  a  matter  of  hours.  A  large  semiannual  effect  in  the  density, 
the  cause  of  which  is  still  not  understood,  was  among  the  variations  modeled.  Some  of 
these  variations  are  illustrated  in  Figs.  2,  3,  and  4. 

In  recent  years  satellites  such  as  Explorer  17,  OV1-15,  Explorer  32,  OGO-6,  ESRO- 
4  and  AE-C,  D,  and  E  have  carried  mass  spectrometers  and  have  therefore  been  able  to 
measure  individual  species.  A  number  of  variations  not  apparent  in  the  satellite  drag 
data  were  discovered  and  included  in  the  most  recent  models,  MSIS  and  Jacchia  1977. 

While  continual  improvements  to  the  models  can  be  anticipated,  there  is  no 
immediate  prospect  of  a  revolutionary  improvement  in  these  models;  thus,  it  is  now  an 
appropriate  time  to  review  modeling  efforts  to  date. 
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SATELLITE  1966  44A  (EXPLORER  52) 
M  JO 

00  800  900  40000  K»  700  XX) 


The  semiannual  density  variation  as  derived  from  drag  analysis  on  Satellite 
1966  44A  (Explorer  32).  All  other  variations  have  been  suppressed  by  using 
the  appropriate  equations  (from  Jacchia,  1971). 


MJO 


Figure  3.  The  geomagnetic  effect  as  derived  from  the  drag  of  six  satellites  in  May 
and  June  1967.  represents  the  change  in  exospheric  temperature 

required  to  represent  the  observed  density.  The  lowest  panel  shows  the 
change  in  exospheric  temperature  computed  from  Kp  by  using  the  model 
(from  COSPAR,  1972). 
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Figure  4.  Atmospheric  temperatures  obtained  on  November  30,  1967,  by  Carru  and 
Waldteufel  (1969)  by  use  of  Thomaon-acatter  techniques,  compared  with 
temperatures  predicted  by  the  J71  models  for  a  height  of  300  km  (solid 
line)  (from  Jacehia,  1971). 


M 


IL  MODEL  DESCRIPTIONS 


Early  Models 

(ARDC  1959;  Jaeehia  1960  and  'DENSEL;  C1RA  1961;  Harris-Priester  and  CIRA  1965; 
Jacchia  1965  and  Jacchia-Walker-tfruee) 

The  basis  and  limitations  of  some  thermospheric  models  that  were  developed  in 
the  decade  following  the  launch  of  the  first  artificial  earth  satellite  are  discussed  in  this 
section.  A  fairly  large  number  of  models  were  developed  in  that  period,  but  the 
discussion  will  be  limited  to  those  that  either  achieved  international  recognition  or  were 
used  frequently  for  Aerospacc/SAMSO  programs.  A  few  of  these  models  are  still  in  use. 
Some  of  the  early  models  attempted  a  "full"  description  of  the  atmosphere;  i.e.  pressure, 
temperature,  and  composition  were  specified  in  addition  to  the  density.  Since,  however, 
such  attempts  were  not  based  on  actual  measurements  of  temperature  or  composition, 
our  discussion  shall  be  limited  only  to  the  density. 


ARDC  1959 

The  ARDC  1959  (Air  Research  and  Defense  Command)  model  (Minzner  et  al., 
1959)  was  the  first  standard  upper-atmospheric  model  that  utilized  satellite  measure¬ 
ments  of  the  density.  Densities  were,  of  course,  deduced  by  the  orbital  decay  method, 
and  so  only  very  coarse  time  and  space  resolution  was  provided  by  the  measurements. 
Accordingly,  only  a  single  profile  of  the  density  was  presented  out  to  700  km  altitude. 
Although  the  data  base  covered  altitudes  from  170  km  to  650  km,  there  was  a  large 


interval  between  about  350  and  600  km  where  there  were  no  measurements.  The 
satellite  orbits  used  were  mostly  of  moderate  inclination  and  covered  different  local 
times;  thus,  the  model  was  representative  of  "average"  global  conditions  Airing  a  period 
of  high  solar  activity  (1957-58).  The  density  profile  above  150  km  was  higher  than  that 
of  the  ARDC  1956  model  (which  was  based  only  on  rocket  data).  Initially  it  was 
suspected  that  the  difference  was  due  to  errors  of  measurement.  Today,  with  our 
increased  knowledge  of  atmospheric  behavior,  we  can  believe  that  the  major  part  of  the 
difference  was  most  probably  related  to  real  solai^eyele  induced  variability  of  the 
atmosphere. 

Jacchia  1960  and  'DENSF.L 

Jacchia*s  (1960)  first  model  of  thermospheric  density  was  based  on  satellite  drag 
data  from  four  satellites  having  perigee  altitudes  between  210  and  650  km.  The  model 
took  into  account  both  diurnal  effects  and  variations  related  to  solar  activity.  Geo¬ 
magnetic  activity  effects  had  been  observed  in  orbital  decay  data  but  were  not  suf¬ 
ficiently  well-understood  to  be  modeled  at  that  time. 

Unlike  the  ARDC  1959  model,  which  was  presented  in  tabular  form,  the  Jacchia 
1960  model  density  was  expressed  analytically  by 


p  =  p0  F20  {  1+0,19  [**P  (°-0055  *  !••]  eo*6  f  } 

m 


(1) 
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where  log  Pq  =  -16.021  -  0.001985  *  ♦  6.363  exp  (-0.0026  z).  In  this  expression,  p  is  the 

3 

density  in  gm/cm  at  altitude  z  in  km  (200<z<700);  F2Q  is  the  20-cm  solar  flux  in  units 
•22  2 

of  10  W/m  -Hz.  represents  the  angular  distance  from  the  maximum  of  the  diurnal 
density  bulge,  which  was  placed  at  the  subsolar  latitude  at  a  local  time  of  about  14  hr. 
(Note  that  Jacchia  used  the  20-cm  solar  flux  as  an  index  of  the  degree  of  solar  activity. 
Subsequent  studies  utilized  the  10.7-cm  flux,  which  was  found  to  be  somewhat  better 
correlated  with  density  variations.) 

The  simplicity  of  the  model  made  it  particularly  useful  for  production-type  cal¬ 
culations  (as  in  orbit  determination)  where  numerical  efficiency  is  required.  The 
Lockheed  Missile  and  Space  Co.  (1961)  recognized  this  and  incorporated  the  Jacchia  1960 
model  with  the  ARDC  1959  model  at  lower  altitudes  for  use  in  orbital  analyses.  The 
Lockheed  product  is  known  as  "Lock heed- Jacchia"  by  users  of  the  Aerospace  TRACE 
program;  it  is  the  same  as  'DENSEL  in  the  Advanced  Orbital  Ephemeris  System  (AOES) 
software. 

In  'DENSEL,  the  density  up  to  an  altitude  of  140  km  is  represented  by  ARDC  1959; 
above  200  km  Jacchia  1960  is  used.  In  the  intervening  region  (140-200  km)  a  smooth 
transition  is  provided  by 


p  =  2.889  x  10 


-12  ^140JL  j  7A8  ' 
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where  p,  z  and  0  have  the  same  meaning  as  in  the  previous  expressions.  It  is  to  be 
noted  that  the  model  utilizes  the  10.7-cm  flux  (Fjq  ^  in  units  of  10~^W/m*-Hz)  rather 
than  the  20-cm  flux.  The  usage  of  the  10.7-cm  flux  is  maintained  even  in  the  region  200- 
700  km,  where  the  Jacchia  1960  formulation  is  adopted.  This  is  done  by  replacing  Fj0  in 
Eq.  (1)  by  0.85  Fj0  where  0.85  is  the  appropriate  scaling  factor  between  the  solar- 
cycle  averages  of  the  20-cm  and  10.7-cm  fluxes.  While  Eq.  (2)  provides  a  smooth 
transition  from  the  ARDC  1959  to  the  Jacchia  1960  model,  the  slope  of  the  density  at 
the  end  points  is  usually  discontinuous.  This  feature  is  generally  not  important  in 
applications. 

At  very  high  orbital  altitudes  (700-1850  km)  'DENSEL  represents  the  density  with 
the  expression 


P  = 


2.60  x  10‘5  F 


10.7 


2  h3 


<» 

h 


where  h  is  the  altitude  in  nmi  (h  =  z/1.852).  Equation  (3)  is  only  approximate,  based  on 
the  very  limited  data  that  were  available  at  the  time.  For  orbital  studies,  however,  it 
generally  is  adequate  since  the  air  <frag  is  very  small  relative  to  other  perturbing  forces. 
In  the  light  of  more  recent  data  we  now  know  that  eq.  (3)  does  not  accurately  describe 
the  global  and  seasonal  variation  of  the  density  (which  is  predominantly  made  up  of  He  in 
that  regime). 
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A  feature  of  'DENSEL  which  should  be  noted  is  that  a  formula  for  the  time- 
dependent  behavior  of  the  10.7-cm  flux  is  used  whenever  actual  values  of  Fjq  ^  are  not 
provided  to  the  program.  This  is  given  by 


Fl0  7  =  150  ♦  80  cos  (2rr  d/4020) 


(4) 


where  d  is  the  time  in  days  from  the  previous  solar-cycle  maximum.  As  can  be  seen, 

-22  2 

Fj0  7  in  Eq.  (4)  varies  between  70  and  230  x  10  W/m  -Hz.  Thus,  the  predicted  values 
encompassing  the  maximum  arc  on  the  high  side.  It  is  to  be  noted  that  'DENSEL  does  not 
include  any  of  the  known  atmospheric  variations  other  than  the  diurnal  and  solar  cycle 
effects. 


CIRA  1961 


The  CIRA  1961  (COSPAR  International  Reference  Atmosphere)  model  consisted  of 
three  tables  of  the  atmospheric  variables  that  represented  the  daily  average  conditions, 
daytime  maximum  and  nighttime  minimum  conditions.  No  other  time  or  space  variation 
of  the  atmosphere  was  described.  The  model  purported  to  be  representative  of  medium 
solar  activity  conditions,  but  since  the  data  base  consisted  of  orbital  decay  observations 
in  the  1958-60  period,  it  should  be  regarded  as  a  high  solar  activity  model. 
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The  lower  boundaries  of  the  three  models  that  comprise  CIRA  1961  are  different. 
The  daily  mean  model  extended  from  ground  level  to  800  km;  the  daytime  maximum 
model  started  at  200  km,  the  nighttime  minimum  model  at  170  km. 

Harris-Priester  and  CIRA  1965 

The  work  of  Harris  and  Priester  (1962a,b)  represents  a  milestone  in  thermospheric 
modeling  studies  since  theirs  was  a  first  major  attempt  to  describe  the  upper  atmosphere 
on  a  physical  basis.  Rather  than  fit  the  observational  data  with  profiles  or  empirical 
formulas  that  had  no  physical  basis,  they  fit  the  data  subject  to  the  constraints  of  the 
thermospheric  energy  equation  under  conditions  of  hydrostatic  equilibrium.  In  other 
words,  they  solved  the  one-dimensional  time-dependent  heat  conduction  equation  to 
obtain  the  temperature,  and  with  this  the  different  constituent  profiles  could  be 
computed  assuming  diffusive  equilibrium.  Once  the  temperature  and  the  concentrations 
of  the  different  species  are  known,  all  other  atmospheric  variables  are  determined.  To 
fit  the  computed  density  to  the  observations,  the  boundary  conditions  or  properties  of  the 
energy  sources  can  be  adjusted. 

Harris  and  Priester  encountered  some  difficulty  in  trying  to  match  the  observa¬ 
tions  using  only  the  solar  EUV  flux  as  the  source  of  heating.  The  amplitude  of  the  daily 
density  variation  was  too  large,  the  time  of  maximum  was  several  hours  too  late,  and  the 
required  heating  efficiency  was  too  high.  All  these  problems  could  be  removed  if  a 
second  heat  source  were  added;  this  hypothetical  source  needed  to  be  similar  to  the  EUV 
source  in  its  strength  and  vertical  distribution,  but  the  maximum  heating  was  required  at 
about  9  a.m.  local  time.  It  was  qualitatively  argued  that  the  second  source  could  be 
related  to  corpuscular  energy  input. 


Using  this  second  heat  source  to  supplement  EUV  heating,  Harris  and  Priester 
(1962a)  produced  a  diurnal  model  of  the  upper  atmosphere  that  was  appropriate  for 
equatorial  and  temperate  latitudes  at  equinox  and  for  solar  activity  conditions 
corresponding  to  ^  =  200.  Shortly  thereafter,  Harris  and  Priester  (1962b)  extended 
their  model  to  include  solar  cycle  variations.  Tables  were  published  in  NASA  document 
D-1444  (1962). 

COSPAR  (Congress  on  Space  Research)  adopted  the  Harris-Priester  models  for  the 
C1RA  1965  model.  Under  this  title  120  tables  were  published  that  described  the 
thermospheric  structure  from  120-800  km  for  10  levels  of  solar  activity  (65  s  FJ0  ^  225) 
and  with  2-hr  increments  of  local  time.  The  tables,  as  described  above,  were  based  on 
theory.  However,  to  account  for  other  observed  density  effects,  such  as  the  2/-day 
variation,  geomagnetic  activity  effects  and  semiannual  changes,  empirically  derived 
corrections  to  the  exospheric  temperature  need  to  be  used  in  conjunction  with  the  tables. 


Jacchia  1965  and  the  Jacchia-Walker-Bruce  Model 

With  satellite  drag  data  from  half  a  solar  cycle  (1958-1964),  Jacchia  (1965)  was 
able  to  construct  an  improved  empirical  model  that  took  into  account  solar  activity 
(including  both  ll-vr  and  27-day  cycles),  diurnal  variations,  and  geomagnetic  activity. 
Although  the  published  model  was  accompanied  by  tables  of  the  atmospheric  quantities, 
the  composition,  temperature  and  density  could  be  directly  computed  for  any  set  of  given 
conditions.  This  is  in  contrast  to  the  CIRA  1965  model,  in  which  one  had  no  choice  but  to 
interpolate  the  tables  for  arbitrary  conditions.  The  formulation  of  the  Jaechia-1965 
model  was  also  continuous  over  the  globe,  whereas  CIRA  1965  was  essentially  only  for 
middle  and  low  latitudes. 
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With  this  model  Jacchia  initiated  the  now  common  practice  of  using  the 
exospheric  temperature  as  a  parameter  to  synthesize  the  density  data  at  different 
altitudes.  To  elaborate,  given  a  lower  boundary  temperature  and  the  exospheric 
temperature,  a  temperature  profile  is  assumed  to  be  determined.  Then,  given  a  set  of 
lower  boundary  conditions  on  number  densities  of  the  constituents,  the  equations  for 
diffusive  equilibrium  are  integrated  to  determine  the  vertical  distributions  of  the 
different  atmospheric  constituents,  which  may  be  summed  to  obtain  the  total  density: 


IfiTz)  <** 


-(1  ♦  a.) 


dT 

T 


(5) 


P  =  ^  njm, 


(6) 


where  n^  is  the  concentration  of  the  ith  constituent,  nij  is  its  mass,  g  is  the  acceleration 
of  gravity,  z  is  the  altitude,  k  is  Boltzmann's  constant,  T  is  the  temperature,  is  the 
thermal  diffusion  coefficient  and  p  is  the  total  mass  density.  In  essence,  then,  given  a 
set  of  lower  boundary  conditions,  a  density  profile  is  determined  by  specifying  the 
exospheric  temperature.  Thus,  rather  than  manipulate  the  density  to  simulate  the 
different  atmospheric  variations,  Jacchia  1965  utilizes  exospheric  temperature  changes 
to  produce  the  density  variations.  The  appropriate  temperature  changes,  which  were 
analytically  formulated,  were  derived  by  fitting  the  model  to  observations. 
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When  Jacchia  published  his  model  in  1965,  it  was  stated  that  the  concentrations 
were  to  be  computed  by  numerical  integration  of  the  diffusion  equation  using  the  appro¬ 
priate  temperature  profile.  Walker  (1965)  pointed  out,  however,  that  with  a  minor 
change  in  the  mathematical  form  of  the  temperature  profile,  the  diffusion  equation  could 
be  integrated  directly.  This  feature  is  quite  important  in  applications  where  numerical 
efficiency  and  economy  are  imperative. 

The  Jacchia-Walker-Bruce  model  (Bruce,  1966),  which  is  a  standard  option  for 
TRACE  users,  is  basically  Jacchia-1965  with  densities  computed  using  the  closed-form 
solution  of  the  integral  of  the  diffusion  equation.  An  additional  modification  introduced 
by  Bruce  eliminates  the  cross-over  of  the  density  profiles  near  150  km  as  shown  in  Fig. 
5.  This  crossing  of  density  profiles,  which  means  that  on  opposite  sides  of  the  cross-over 
(or  isopynic)  point  the  density  varies  in  opposite  senses  as  a  function  of  exospheric 
temperature,  is  due  to  the  conservation  of  mass.  In  other  words,  the  density  at  high 
altitudes  can  increase  only  at  the  expense  of  the  density  at  low  altitudes  if  the  pressure 
at  the  lower  boundary  is  to  remain  constant.  While  it  is  quite  likely  that  such  behavior 
does  exist  somewhere  in  the  atmosphere,  observations  indicate  that  it  is  below  about  the 
120-km  level.  Thus,  Bruce's  modification,  while  leading  to  nonconservation  of  mass,  or 
equivalently,  flow  across  the  boundary,  does  force  the  model  to  behave  in  a  more 


realistic  fashion. 


CIRA  1972/JACCH1A  1971 


Introduction.  The  COSPAR  International  Reference  Atmosphere  1972,  or  CIRA 
1972,  (COSPAR,  1972)  consists  principally  of  two  independent  models;  A  model  of  G.  V. 
Groves  for  the  altitude  range  25-125  km  and  a  model  of  L.  G.  Jacchia  for  the  range  110- 
2000  km.  The  two  models  differ  considerably,  including  the  use  of  different  parameters; 
they  also  do  not  match  in  the  region  of  overlap.  Our  concern  in  this  report  is  with  the 
high  altitude  model  of  Jacchia,  which  is  essentially  identical  with  the  model  commonly 
referred  to  as  Jacchia  1971  or  J71  (Jacchia,  1971). 

The  model  is  based  primarily  on  satellite-drag  data,  with  composition  adjusted  to 
reflect  theoretical  results  and  in  situ  mass  spectrometer  results.  In  particular  the 
controversial  analysis  of  atomic  oxygen  measurements  by  von  Zahn  (1970)  has  been 
incorporated.  As  a  static  model  which  assumes  diffusive  equilibrium,  J71  builds  on  the 
foundations  laid  by  Nicolct  (1965),  Jacchia  (1965),  and  Jacchia  (1970).  Perhaps  the  most 
notable  change  in  J71  is  in  moving  the  constant  boundary  conditions  from  120  km  down  to 
90  km,  which  has  the  effect  of  moving  the  resultant  notorious  isopycnic  layer  from  near 
150  km  down  to  a  height  well  below  normal  satellite  altitudes,  thus  making  the  Bruce- 
type  correction  in  a  more  self-consistent  fashion. 

The  150  km  level  is  employed  in  effect  as  a  reference  level;  the  temperature 
profile  at  lower  altitudes  is  "adjusted  to  give  the  required  total  density  at  150  km" 
(COSPAR,  1972)  after  setting  the  composition  at  150  to  correspond  to  the  composition 
suggested  by  von  Zahn  (1970). 

The  CIRA  1972  model  also  represents  a  model  different  from  the  more  "mathe¬ 
matical"  models,  such  as  MSIS  (to  be  discussed  below),  for  which  the  choice  of  functional 
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forms  is  dictated  by  questions  of  orthogonality,  completeness,  etc.,  with  the  final  para¬ 
meters  determined  by  a  fitting  routine.  In  the  CIRA  1972  mcdel  the  mathematical  forms 
and  the  specific  parameters  enter  in  ways  which  attempt  to  reflect  the  current 
understanding  of  the  physics  and  chemistry  of  the  atmosphere,  or  which  specifically 
demonstrates  the  form  of  a  given  variation,  thereby  isolating  certain  features  or 
variables. 

Model  Formulation.  The  J71  model  is  a  static  model  which  assumes  diffusive 
equilibrium.  The  basic  independent  variable  is  the  temperature.  All  temperature  pro¬ 
files  begin  from  a  constant  value  T^  =  183°K  at  the  base  altitude  of  90  km  and  have  an 
inflection  point  at  125  km.  The  choice  of  the  functional  form  of  the  profile  is  dictated 
as  much  by  computational  as  by  theoretical  and  experimental  considerations.  The  chosen 
profiles  yield  continuous  first  and  second  derivatives  with  respect  to  height  through  the 
inflection  point.  The  temperature-height  profile  is  characterized  by  the  exospheric 
temperature,  Tx  ,  the  temperature  which  is  approached  asymptotically  at  high  altitudes. 

The  J71  model  recognizes  several  types  of  variations  which  are  dealt  with 
explicitly.  Specifically  they  are: 

1.  Solar  cycle  variaions  and  shorter-term  solar  activity  effects. 

2.  The  diurnal  variation. 

3.  Geomagnetic  activity  variations. 

4.  The  semiannual  variation. 

5.  Seasonal-latitudinal  variations. 


(More  rapid  density  variations  (waves)  are  also  recognized  to  be  present  but  as  being 
incapable  of  being  handled  in  a  static  model  of  this  type.) 

The  above  five  sources  of  variations  are  accounted  for  in  the  model  in  two  distinct 
ways.  The  first  is  as  a  variation  of  the  exospheric  temperature  T,,  .  The  variations  so 
modeled  are: 

1.  The  solar  activity  effects,  which  employ  for  their  input  parameters  the 
readily  available  10.7-cm  solar  flux  (FjQ  7).  The  model  distinguishes 
between  a  contribution  from  active  areas  and  a  disk  component,  which  are 
modeled  by  a  slightly  delayed  (1.71  days)  value  of  F10.7  and  by  a  running 
average  over  6  solar  rotations  (Fjfl  7),  respectively. 

2.  The  diurnal  variation  is  represented  by  a  sub-solar  "bulge."  The  ratio  of  the 
maximum  exospheric  temperature  at  the  sub-solar  {joint  to  the  minimum  at 
the  anti-solar  point  is  a  constant  1.3,  although  Jacchia  points  out  the 
considerable  uncertainty  on  this  point.  The  diurnal  variation  has  both  local 
time  and  latitude  dependence  (and  seasonal  variability  implicitly  through 
the  movement  of  the  sub-solar  point).  The  hours  of  minimum  and  maximum 
of  the  daily  density  are  independent  of  latitude  and  altitude  and  occur  at 
02.87  and  14.08  hours  local  solar  time,  respectively. 

The  second  way  in  which  the  sources  of  variation  are  incorporated  into  the  model 
is  via  a  direct  percentage  variation  of  the  local  ambient  density.  The  sources  modeled  in 
this  way  are: 


1. 


The  semiannual  variation  is  expressed  as  a  separable  function  of  time 
(expressed  in  days)  and  altitude  in  the  manner 


"*10  ^semiannual  *  (7) 


The  function  g(t),  which  is  not  sinusoidal,  has  maximum  excursions  of 
-0.519  and  *0.478  and  averages  to  zero  over  the  period  of  a  year.  The 
height  function  ftz)  is  relatively  small  below  200  km  with  a  broad  maximum 
peaking  at  0.353  near  800  km. 

2.  The  seasonal-latitudinal  variations  in  the  lower  thermosphere  arise  in  J71 
as  an  attempt  to  fit  lower  boundary  conditions.  The  form  chosen  is 

4  '°«10  f>  *  8  -f-  P  Si"2*  (1) 


where  0  is  the  geographic  latitude,  S  is  a  function  of  altitude  which  has  a 
maximum  value  of  0.166  at  110  km,  falling  to  zero  by  170  km,  and  P  is  a 
sinusoidal  function  with  the  maximum  at  the  winter  solstice. 

Another  seasonal-latitudinal  variation,  namely  that  of  helium,  is  modeled 
separately  as  a  function  of  solar  declination  and  latitude.  The  maximum, 
d  log]Q  n  (He)  =  0.420,  occurs  over  the  winter  pole  at  solstice.  The  summer  polar 
minimum  is  -0.230. 
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The  variations  with  geomagnetic  activity  are  modeled  by  a  combination  of 


temperature  and  density  variations.  The  suggested  formulas  are: 

4l°gio  P=  0.012  Kp  ♦  1.2  x  10'5  exp(Kp) 


(9) 


AT  *  14°K  ♦  rf?02exp(K  >. 

<s  p  P 

An  alternative  version  employs  only  a  temperature  variation 

AT a  =  28°Kp  ♦  (f?03  exp  (Kp)  . 


(10) 


(11) 


A  time  lag  of  6.7  hours  for  the  appropriate  value  of  Kp  is  suggested. 

In  practice,  the  model  is  used  by  fit  it  determining  the  appropriate  jxospheric 
temperature,  T^  ,  according  to  the  jrious  parameters  (local  time,  time  of  year,  Fj0 
etc.)  which  directly  affect  this  quantity.  The  value  of  T^  then  completely  determines  a 
given  density-height  profile  for  the  individual  species  according  to  the  equations  for 
diffusive  equilibrium.  Adjustments  to  these  individual  density  profiles  are  then  made 
according  to  the  formulas  for  the  semiannual,  geomagnetic,  and  seasonal-latitudinal 
effects  which  operate  directly  on  the  density. 

Limitations.  CIRA  1972/J71  is  a  static  model  derived  primarily  from  insights 
gained  from  examination  of  satellite  drag  data;  its  limitations  derive  in  large  part  from 
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these  circumstances.  Because  it  is  a  static  model  it  is  best  able  to  represent 
atmospheric  variability  with  relatively  long  time  scales  like,  for  example,  the  semiannual 
variation.  Variations  with  shorter  time  scales,  such  as  the  diurnal  variation,  are  modeled 
less  well.  For  example,  the  diurnal  response  of  various  constituents  in  the  upper 
atmosphere,  in  particular  Nj  and  O,  have  different  phases  resulting  in  a  change  in  phase 
of  the  overall  diurnal  response  as  the  atmosphere  changes  from  a  dominantly  to  a 
dominantly  O  atmosphere  with  increasing  height.  While  J71  does  contain  diurnal 
asymmetry,  it  has  a  very  simple  form,  constant  at  all  altitudes,  and  is  subject  to 
considerable  uncertainty.  At  even  shorter  time  scales  the  more  dynamic  response  is 
still  more  poorly  modeled,  particularly  for  high  values.  This  difficulty  is  exacerbated 
when  the  comparison  is  made  with  in  situ  mass  or  density  measurements  with  their 
inherently  more  rapid  response  than  the  drag  data  used  for  the  model  formulation. 

Accuracy.  Figure  fi, taken  from  the  CIRA  1972  publication  (COSPAR,  1972),  gives 
density  residuals  from  the  model  for  nine  satellites.  The  residuals  appear  to  average 
about  15-25%  over  the  time  period  1958-1970,  i.e.,  an  entire  solar  cycle.  This  is  in 
agreement  with  values  derived  comparing  the  J71  model  with  in  situ  ion  gauge  data 
below  300  km  for  the  period  1973-1975  (Rice  and  Sharp,  unpublished).  Because  of  lack  of 
available  low  altitude  data  when  formulating  J71,  the  model  at  low  altitudes  may  not  be 
correct  in  all  its  details.  For  example,  the  phase  of  the  diurnal  variation  changes  at  low 
altitude,  a  feature  which  is  not  modeled  by  J71  (Sharp  et  al.,  1978).  Because  the 
amplitude  of  such  variations  also  decreases  at  low  altitudes,  however,  the  effect  on  the 
overall  percentage  accuracy  is  not  extremely  great.  Nevertheless,  care  must  be 
exercised  at  low  altitudes  if  specific  details  such  as  diurnal  behavior  or  relative 
composition  are  significant  for  the  use  considered. 


-30- 


*J0 


ICMO  2  I  H4  041 

I  rr~J 

_ — 

j 

i/r  v  »  • 

•  oi  r_ 

?  •  TOO  -  ©*>  M> 

-Of  l  - ± - 


•  Of 

-..V  ' 

.  V. 

—  •• 

0 

•Of 

r^-73 

t  _  _  __  1 

•  Of 

>**GuAao  7 

1  a  ||  ...  ^  * 

w 

•Of 

I'vJ*  ■  '1  » 

1  •  4*4  w»  ■  .  * 

r.  _  __  .1 

•J  •  #*  •  •* .  '**•  \ 


wtNtuMO  i  irv*  #  >) 
I  .  T>0  »">  i 


O 

r 

<3  *0» 

0 

•  Of  l 


1  HI'W  IW  _»*■  /**■**•.•' 

I  •  MO-  MO  *m  *r{.  *  • «  y 


•  Of 
0 

•Of 


•  Of 

0 

•0* 

•  0» 

0 

•0 
•  Oi 


C«PU*f»  • 


,  |.-***«-*l  \  \  •  1_  T  |.>T»*| 

-WSN-x-'-  j'_  v-  j 

>( — pawwaijiiiMM  T~iqo»» -  {  ■.■■.•,•*'• J 

IT  wio**  — VV.*~-  *1 

K  SST*  T  *u  T  *•’  T  *m  ’  <"*T  •’o 


0)  tHj 
O'  1  ~T 


«*A0I  5  JfOM  fli 


Figure  6.  Thirty-day  means  of  density  residuals  from  the  J71/CIRA  1972  model  for 
nine  satellites  covering  a  full  solar  cycle,  r.  is  an  average  effective  height 
for  each  satellite  (from  COSPAR,  1972). 
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The  McDonnell  Douglas  Model,  MDAC 

MDAC  is  an  atmospheric  model  developed  by  W.  P.  Olson  et  al.  (1975).  The 
motivation  for  developing  the  MDAC  model  arose  as  a  result  of  studies  of  the  dayside 
cusps  which  indicated  that  at  high  latitudes  the  thermospheric  densities  arc  quite 
variable,  with  particle  precipitation  being  the  dominant  energy  source.  Specifically  this 
precipitation,  which  is  controlled  by  the  magnetospheric  magnetic  field,  consists 
principally  of  solar  wind  protons  and  electrons  which  precipitate  directly  into  the 
ionosphere  and  upper  atmosphere  near  noon  (Frank,  1971;  Heikkila,  1971).  No  existing 
thermospheric  model  attempts  to  incorporate  this  particular  energy  source. 

Prior  to  construction  of  the  MDAC  model,  extensive  work  was  done  on  energy 
sources  for  the  thermosphere  (Moe,  1971;  Olson,  1971).  In  addition,  the  physics  and 
morphology  of  the  day-side  cusp  were  investigated  using  available  data  on  precipitating 
particles.  The  intersection  of  the  cusp  is  assumed  to  be  at  75°  magnetic  latitude  with  its 
central  longitude  on  the  magnetic  noon  meridian  plane.  The  extent  of  the  region  of 
interaction  Is  assumed  to  be  180°  in  longitude  and  3-5°  in  latitude. 

The  MDAC  model  is  semi-empirical  and  has  a  range  of  0-450  km.  Between  120 
and  450  km  the  model  provides  a  global  description  of  the  atmosphere,  while  below 
120  km  it  utilizes  a  power  series  fit  to  the  CIRA  1972  mean  reference  atmosphere.  The 
boundary  between  the  lower-  and  upper-atmospheric  models  is  everywhere  continuous 
and  the  slopes  are  closely  matched. 

The  total  density  returned  by  the  model  is  the  sum  of  two  terms:  the  corpuscular 
heating  term  (  A p^)  is  described  in  magnetic  coordinates,  while  the  EUV  and  other 
sources  term  (  )  is  described  in  geographical  coordinates.  The  two  contributions  are 
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then  •'nmbined  in  the  following  manner: 

P  =  P„  *  APC  -  P(*.  A,  ^  ,  T,  D,  F,  P)  in  gm/em3 

z  =  altitude  in  km 

A  =  geographic  latitude  in  degrees 

0  =  east  geographic  longitude  in  degrees 

T  =  universal  time  in  hours 

D  =  day  of  the  year 

-22  2 

F  =  deci metric  solar  flux  in  10  W/m  -Hz 

P  =  particle  heating  parameters. 

The  first  term  is  expressed  as  a  product 

P„  =  p0  (z,  F)  B  (z,  X  ,  t)  J  (z,  X  ,  O)  Q  (z,  D) 

where 

PQ  is  the  mean  equatorial  density 

B  is  the  diurnal  variation 
J  is  the  seasonal-latitudinal  variation 


Q  is  the  semi-annual  variation 


t  is  the  local  time. 


For  average  conditions  pQ  can  be  expressed  as 


Po  =  p(120)  exp  - 120  -  where  p(120)  =  2.7  x  10"11  gm/cm3.  (14) 

1.737 v/  z-103 


The  F10  7  variation  is  introduced  by  replacing  1.737  with 

0.99  ♦  0.518(  * 


where  F  is  averaged  over  the  preceding  3  months.  The  diurnal  variation  is  given  by 


B(z,  X  ,  t)  =  1  ♦  (f(t)  -  1)  cos  X 


The  function  ft t)  is  the  diurnal  variation  in  density  at  the  equator  given  by 


f(t)  =  0.994  ♦  0.545  eos^—  (t  -  14.745)j  ♦  0.102  cos  |^-  (t  -  1.838)J  .  (17) 


The  function  ^(^)  is  given  by  1.1  x  (1  -  exp(-(z-120)/150))  which  satisfies  the  following 
conditions.  It  reproduces  the  values  at  450  km  observed  by  Hedin  et  al.  (1974).  The 
height  dependence  fits  the  LOGACS  data  (Bruce,  1972)  near  local  noon  at  188  km  and 
the  data  of  Ching  and  Carter  (1974)  near  local  midnight  at  230  km.  In  addition  the 
diurnal  variation  is  made  to  approach  zero  where  it  interfaces  with  CIRA  1972  at 


120  km. 


The  latitudinal-seasonal  dependence  is  given  by 

m  /i(a) 

J  (z,  X,D)  =  1  -  a(l  -  eos(  X-  Xs» 

Xs  is  the  latitude  of  the  subsolar  point,  23.5  sin  deg. 

a  =  0.25  and  m  =  2.5. 

This  is  a  compromise  equation  to  fit  the  latitude  dependence  observed  by  Hedin  et  al. 
(1974)  at  450  km,  the  LOGACS  data  at  186  km,  and  observations  indicating  only  a  small 
latitude  dependence  due  to  UV  heating. 

The  semiannual  variation  is  expressed  as 

Q(z,  D)  =  1  ♦  R  (z)  G  (D)  (20) 


(18) 

(19) 
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The  corpuscular  heating  term  is  given  by 


A  Pe  =  e*P  [-  (z-120)/50]c  (  Am,  tm,  P) 


(23) 


where  Am  is  the  magnetic  latitude,  t  is  the  magnetic  local  time  and  P  the  cusp  heating 
parameters  which  include  the  particle  precipitation  intensity  and  size  parameters  of  the 
heating  region.  The  function  is  given  a  scale  height  of  50  km  in  order  to  have  zero 
heating  at  120  km  and  approach  a  maximum  heating  near  170  km.  This  is  done  because, 
although  protons  can  penetrate  to  about  110  km,  the  peak  energy  loss  per  km  is  at 
-169  km  (Olson,  1971,  1972). 

Function  C  represents  a  density  bulge  which  encircles  the  geomagnetic  pole  and  is 
most  intense  at  local  magnetic  noon.  The  amplitude,  extent  and  location  of  the  bulge  are 
adjustable  parameters  which  depend  on  solar  wind  and  magnctospheric  characteristics. 
These  adjustable  values  have  been  nominally  set  to  give  a  best  fit  to  the  LOOJACS  data 
during  quiet  times. 

The  MDAC  model  has  the  advantage  of  being  functionally  simple  and  compu¬ 
tationally  fast.  The  FORTRAN  version  of  this  model  (as  provided  by  Olson,  et  al.)  has 
been  carefully  optimized  to  minimize  the  number  of  multiplications  and  calls  to 
elementary  functions.  Since  the  model  does  not  include  the  Kp  variation,  its  utility  is 
limited  to  quiet  times.  In  addition  it  predicts  densities  which  are  in  general  about  10% 
too  low.  However,  since  it  does  include  the  heat  source  from  charged  particle 
precipitation  into  the  high  latitude  regions,  it  can  be  valuable  as  a  baseline  in  the  study 
and  prediction  of  densities  in  the  polar  thermosphere. 


-36- 


The  U.  S.  Standard  Atmosphere,  1976 

Introduction.  The  U.  S.  Standard  Atmosphere,  1976,  is  the  atmospheric  model 
adopted  officially  by  the  United  States  Committee  on  Extension  to  the  Standard  Atmo¬ 
sphere  (COESA).  It  replaces  an  earlier  model  known  as  "U.S.  Standard  Atmosphere, 
1962"  (COESA,  1962).  The  new  standard  is  defined  and  described  with  graphs  and  com¬ 
plete  tables  in  the  official  publication  "U.  S.  Standard  Atmosphere,  1976"  (COESA,  1976) 
published  jointly  by  the  National  Oceanic  and  Atmospheric  Administration  (NOAA),  the 
National  Aeronautics  and  Space  Administration  (NASA),  and  the  United  States  Air  Force. 
The  new  standard  is  also  described  together  with  abbreviated  tables  by  Minzner  et  al. 
(1976). 

The  U.  S.  Standard  Atmosphere,  1976,  attempts  to  set  forth  an  idealized  repre- 

|  • 

sentation  of  the  atmosphere  which  will  have  utility  for  a  number  of  scientific, 
engineering,  meteorological,  and  commercial  needs.  As  such  it  differs  from  the  other 
recent  models  discussed  in  this  document  in  two  significant  ways.  Firstly,  only  a  single 
model  profile  is  produced  which  is  "designed  to  represent  mean  mid-latitude  conditions  as 
presently  known  for  heights  within  and  below  the  thermosphere  and  to  represent  mean 
solar  conditions  in  and  below  the  exosphere"  (COESA,  1976).  Secondly,  the  altitude  range 
covered  by  the  model  extends  from  ground  level  (actually  from  5000  meters  below  sea 
level)  upwards  to  exospheric  altitudes  (1000  km).  The  complete  CIRA  1965  and  CIRA 
1972  models  do  not  extend  below  30  km  and  25  km,  respectively. 

Both  of  the  above  unique  characteristics  derive  from  the  manner  in  which  the 
concept  of  a  "standard  atmosphere”  has  evolved.  The  century-long  history  of  the 
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development  of  and  changes  in  "standard"  model  atmosphere  is  well  told  by  Minzner 
(1977).  The  initial  impetus  for  their  use  was  as  a  calibration  standard  for  aneroid 
barometers  used  as  height  measuring  instruments.  The  subsequent  development  of 
aviation  and  its  rapid  expansion  following  World  War  I  increased  the  need  for  and  the 
accuracy  required  of  such  models.  Development  of  higher  flying  aircraft  and  rockets 
during  and  following  World  War  II  greatly  increased  the  «ed  for  a  higher  altitude  model 
and  at  the  same  time  increased  the  data  base  available  for  its  construction.  The 
launching  of  Sputnik  I  in  1957  and  the  subsequent  proliferation  of  satellites  circling  the 
earth  produced  a  similar  situation  extending  out  to  the  exosphere.  The  U.  S.  Standard 
Atmosphere,  1962,  incorporated  much  of  this  early  satellite  data  in  producing  a  model 
which  reached  to  700  km  altitude.  The  primary  impetus  for  issuing  the  U.  S.  Standard 
Atmosphere,  1976,  was  recognition  that  the  "mean"  conditions  chosen  for  the  1962  model 
were,  in  fact,  strongly  biased  toward  conditions  representative  of  high  solar  activity.  As 
a  result,  the  thermospheric  temperature  profile,  in  particular,  was  considerably  hotter 
and  hence  densities  considerably  greater  than  subsequent  measurements  have  supported. 
Differences  between  the  1962  and  1976  standards  are  confined  to  altitudes  above  51  km. 
In  this  report,  our  interest  is  with  the  upper  atmosphere,  which  for  this  model  can  be 
defined  as  the  region  above  86  km  where  a  change  in  computational  techniques  is  made. 
The  following  discussion  will  refer  only  to  this  higher  region  unless  stated  otherwise. 

Model  Formulation.  The  U.  S.  Standard  Atmosphere,  1976,  is  intended  to  repre¬ 
sent  year-round  mean  atmospheric  properties  for  a  mid-latitude  location  during  moderate 
solar  conditions.  Quoting  the  COESA  Working  Group  with  respect  to  the  upper  altitudes: 
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"The  atmosphere  shall  also  be  considered  to  rotate  with  the  earth,  and  be  an  average 
over  the  diurnal  cycle,  semi-annual  variation,  and  the  range  of  conditions  from  active  to 
quiet  geomagnetic,  and  active  to  quiet  sunspot  conditions.  Above  the  turbopause  (above 
110  km)  generalized  forms  of  the  hydrostatic  equations  apply"  (COESA,  1976).  This 
definition  leaves  open  the  determination  of  just  what  average  conditions  actually  obtain 
in  the  earth's  atmosphere,  a  difficulty  which  grows  with  increasing  altitude.  This  subject 
will  be  discussed  further  when  we  examine  the  model's  limitations. 

The  primary  parameter  for  defining  the  U.  S.  Standard  Atmosphere,  1976,  is  the 
temperature  vs.  altitude  profile,  as  it  has  been  for  every  accepted  reference  atmosphere 
since  the  constant  temperature  (50°F)  British  standard  atmosphere  published  by  the 
Astronomer  Royal  G.  B.  Airy  in  1867.  The  temperature  profile  is  used  to  achieve 
continuity  between  the  various  altitude  regimes  where  differing  computational  tech¬ 
niques  must  be  employed.  The  actual  form  of  the  temperature  profile  chosen  must  be  a 
reasonably  accurate  representation  of  the  average  measured  profile  and  at  the  s  me  time 
be  of  a  form  suitable  for  use  with  the  computational  formulas.  Prior  to  the  1976  model, 
the  temperature  profiles  consisted  of  segments  linear  in  either  temperature  T  or 
molecular  scale  temperature  MT  with  respect  to  either  geometric  or  geopotential  height. 
The  1976  model  also  employs  linear  temperature  segments  up  to  86  km.  Above  86  km  the 
1976  standard  uses  a  temperature  profile  consisting  of  four  (non-linear)  segments  chosen 
so  that  the  first  derivative  of  temperature  with  respect  to  geometric  height  is  continuous 
over  the  entire  region  86-1000  km.  The  exospheric  temperature  which  is  approached 
asymptotically  at  high  altitude  is  1000°K,  500°K  below  the  exospheric  temperature  used 
in  the  1962  standard.  The  1962  and  1976  profiles  of  kinetic  temperature  are  shown  in 
Figure  7. 
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Figure  7. 


Kinetic  temperature  vs  geometric  altitude  for  the  1962  and  1976  Standard 
Atmospheres.  Note  in  particular  the  differences  near  100  km  and  above 
120  km  (from  CO  ESA,  1976). 


The  computation"!  are  performed  differently  in  the  regions  below  and  above  86  km. 
Below  the  86  km  the  model  assumes  ground  level  mixing  ratios  (except  for  trace  con¬ 
stituents  such  as  ozone  which  are  handled  separately)  and  hydrostatic  equilbrium.  Above 
86  km  the  simplified  hydrostatic  equation  no  longer  applies,  and  diffusion  and  vertical 
transport  as  well  as  effects  of  atmospheric  chemistry  mast  be  considered.  Thr  transition 
to  a  purely  di.fusive  regime  is  accomplished  gradually  in  the  region  86-115  km.  It  is 
therefore  necessary  to  consider  each  constituent  separately;  many  of  the  atmospheric 
properties  which  are  calculated  directly  below  86  km,  such  as  pressure,  viscosity,  etc., 
mast  be  determined  as  the  sum  of  the  contributions  of  the  various  constituents  at  higher 
altitudes.  A  number  of  processes  known  to  occur  in  this  height  range  are  taken  into 
account  by  a  combination  of  tractable  analytical  formulas  and  empirically  determined 
constants  to  reproduce  experimental  data.  Since  small  adjustments  in  parameters 
propagate  their  effects  to  higher  altitudes,  the  choice  of  a  particular  parameter  value  (or 
even  functional  form)  is  often  the  result  of  an  iterative  procedure.  The  processes 
considered  in  formulating  the  model  are 

1.  Vertical  transport 

2.  Height-dependent  molecular  diffasion 

3.  Height-dependent  eddy  diffusion 

4.  Thermal  diffusion 

5.  Dissociation  and  recombination. 
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The  altitude  profiles  for  the  various  species  considered,  N^,  O,  02<  A r,  He,  and  H, 
should  be  determined  simultaneously  as  a  coupled  set  of  equations  using  the  predeter¬ 
mined  temperature  profile.  In  practice,  in  order  to  reduce  computation  time,  an 
iterative  procedure  was  used,  starting  with  and  proceeding  in  the  order  in  which  the 
species  are  named  above.  The  rationale  for  this  is  explained  in  Part  1  of  the  published 
model  (COESA,  1976).  The  parameters  chosen  for  the  1976  model  have  the  eddy 
diffusion  coefficient  set  to  zero  above  115  km,  the  thermal  diffusion  and  vertical 
velocity  terms  small  or  zero  above  120  km,  and  complete  diffusive  equilibrium  obtaining 
above  150  km.  Hydrogen  is  an  exception  to  the  above;  it  is  only  determined  for  altitudes 
150-1000  km,  with  diffusive  equilibrium  obtaining  essentially  above  500  km. 

After  the  number-density  altitude  profiles  have  been  obtained  for  the  individual 
species,  a  number  of  other  atmospheric  properties  are  determined  by  combining  the 
profiles  with  the  previously  determined  temperature  profile.  These  include  pressure, 
total  number  and  mass  densities,  density  and  pressure  scale  heights,  mean  air-particle 
speed,  mean  free  path,  mean  collision  frequency,  speed  of  sound,  dynamic  and  kinematic 
viscosity,  and  the  coefficient  of  thermal  conductivity.  All  of  these  parameters  are 
either  tabulated  or  shown  graphically  in  the  published  model.  The  density-height  profile 
for  the  1976  standard  and  several  earlier  standards  are  shown  in  Figure  8. 

Accuracy  and  Limitations.  When  considering  a  reference  atmosphere  such  as  the 
1976  U.S.  Standard  the  question  of  accuracy  involves  somewhat  different  considers tions 
than  for  a  more  complete  model  such  as  CIRA  1972  or  Jaechia  1977. 
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Density-height  profiles  of  several  model  and  standard  atmospheres  showing 
the  evolution  from  1947  to  1976.  Note  in  particular  altitudes  above  about 
250  km  where  the  range  is  more  than  an  order  of  magnitude  (from  Minzner, 


a.  Appropriate  Mean  Conditions 

Since  the  model  is  meant  to  represent  mean  conditions,  the  first  question  must  be: 
How  appropriately  are  the  mean  conditions  selected  and  how  well  do  the  selected  data 
represent  this  mean  condition.  .As  noted  previously,  a  primary  impetus  for  issuing  the 
1976  Standard  was  the  recognition  that  the  "average"  conditions  chosen  for  the  1962 
U.  S.  Standard  Atmosphere  were  strongly  biased  toward  conditions  of  high  solar  activity. 
The  1976  Standard  is  based  on  satellite  data  covering  almost  two  complete  solar  cycles 
and  rocket  data  covering  almost  three  complete  solar  cycles;  it  is  therefore  probable 
that  the  1000°K  exospheric  temperature  and  details  of  the  profile  fairly  closely 
approximate  true  "mean"  conditions.  Nevertheless,  future  data  may  force  re-evaluation. 
The  problem  is  also  complicated  by  the  range  of  variability  found  in  the  atmosphere.  At 
lower  altitudes  this  range  is  not  great,  but  at  higher  altitudes  the  deviations  increase 
dramatically,  as  can  be  seen  in  the  temperature  profiles  in  Figure  9.  With  such  large 
variations,  the  choice  of  an  appropriate  "mean"  profile  becomes  as  much  a  matter  of 
taste  or  utility  as  of  science. 

b.  Data  Base  Accuracy 

Considering  the  great  variability  in  the  upper  atmosphere,  it  may  be  that 
questions  regarding  the  accuracy  of  the  data  used  to  construct  the  1976  Standard  are  not 
crucial.  Nevertheless  such  questions  are  particularly  important  in  the  range  86-150  km 
where  the  transition  to  a  diffusive  model  takes  place.  Much  of  the  chemistry  and  the 
dynamics  of  this  region  is  not  well-known  and  in  many  cases  the  few  available 
measurements  appear  to  conflict  with  requirements  imposed  by  measurements  at  higher 
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u*10„f<^pst(t) 

Figure  9.  Range  of  deviations  of  the  density-altitude  profiles  from  the  1976  standard 
for  various  degrees  of  solar  activity:  (A)  for  lowest  exospheric  tempera¬ 
ture  at  sunspot  minimum;  (B)  average  conditions  at  sunspot  minimum;  (C) 
average  conditions  at  average  sunspot  maximum;  and  (D)  exceptionally 
active  conditions  (from  COESA,  1976). 
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and  lower  altitudes  and  current  theoretical  understanding.  Future  changes  to  the 
Standard  will  probably  have  their  origin  in  improved  understanding  of  this  critical  region, 
with  effects  possibly  propagating  in  the  model  to  higher  altitudes  as  well. 

c.  Scientific  Validity 

Anyone  contemplating  use  of  the  1976  U.  S.  Standard  Atmosphere  for  compu¬ 
tational  or  for  theoretical  purposes  must  consider  how  well  the  model  profile 
approximates  a  valid,  physically  possible  profile,  both  for  the  individual  constituents  and 
as  an  ensemble.  COESA  has  attempted  to  facilitate  such  use  by  selecting  differentiable 
temperature  profiles.  As  before,  the  crucial  region  is  the  least  understood,  namely  85- 
150  km  altitude.  Even  a  mean  atmosphere  should  ideally  represent  a  condition  which  can 
actually  obtain  in  the  real  atmosphere.  Because  of  the  highly  dynamic  nature  of  the 
atmosphere,  it  is  possible  that  no  static  model  can  fulfill  this  criterion  completely. 

In  summary,  it  is  the  utility  and  usefulness  of  the  U.  S.  Standard  Atmosphere  1976 
and  not  any  artificially  defined  "accuracy"  which  mast  determine  its  applicability  in  a 
given  instance.  It  is  useful  because  it  is  a  "standard"  and  therefore  uniquely  determined, 
but  in  a  highly  varying,  highly  dynamic  atmosphere  opportunities  for  its  use  are  limited. 
When  issuing  the  1962  U.  S.  Standard  Atmosphere,  COESA,  with  admirable  candor, 
designated  the  height  region  32-90  km  "tentative"  and  the  region  90-700  km 
"speculative."  Even  15  years  later  it  is  well  to  remember  the  warning  implicit  in  these 
words. 
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The  OGO-6  and  ESRO-4  Models 


Introduction.  Two  models  of  thermospheric  temperature,  density  and  composition 
have  been  derived  entirely  from  data  taken  by  satellite-borne  mass  spectrometers.  Since 
these  models  are  identical  in  their  mathematical  formulation,  differing  only  in  their  data 
bases,  they  will  be  described  together. 

Hedin  et  al.  (1974)  developed  a  global  model  of  the  thermosphere  above  120  km 
for  quiet  geomagnetic  conditions  from  measurements  of  the  density  of  Ng,  O  and  He 
taken  by  the  mass  spectrometer  on  the  OGO-6  satellite.  This  satellite  was  launched  on  5 
June  1969  into  an  82°  inclination  orbit  with  398  km  perigee  and  1100  km  apogee  and 
operated  until  July  1971.  During  this  period,  the  latitude  of  perigee  of  the  satellite 
moved  six  times  around  the  earth  and  through  ten  diurnal  cycles.  Values  of  the  mean 
Fj0  7  flux  varied  from  108  to  168,  the  overall  average  being  near  150;  the  average  value 

of  A  was  4.  The  data  used  in  the  model  covered  the  altitude  range  400-600  km. 

P 

Data  taken  by  the  mass  spectrometer  on  board  the  ESRO-4  satellite  during  geo- 
magnetically  quiet  times  were  used  by  von  Zahn  et  al.  (1977)  to  develop  a  global  model  of 
the  thermosphere  above  120  km.  The  ESRO-4  satellite  was  launched  on  22  November 
1972  into  a  91°  inclination  orbit  with  initial  perigee  at  245  km  and  apogee  at  1177  km 
and  operated  until  15  April  1974.  Data  from  228  days  during  this  period  were  used  as  the 
data  base.  Values  of  F1Q  7  ranged  from  ~70  to  M25,  with  a  mean  value  of  92;  the 
average  value  of  Ap  was  7.  The  data  used  in  the  model  covered  the  altitude  range  240- 
320  km.  The  ESRO-4  model  thus  complements  the  OGO-6  model  with  measurements  at 
lower  altitudes  and  lower  levels  of  solar  activity. 
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Model  Formulation.  In  both  models,  the  composition  is  described  in  terms  of 
density  variations  at  a  reference  altitude  zQ  (450  km  for  OGO-6  and  275  km  for  ESRO-4) 
with  an  altitude  dependence  which  assumes  isothermal  conditions  at  higher  altitudes  at 
any  given  time  and  location.  A  spherical  harmonic  expansion  (in  geographic  latitude- 
local  time  coordinates)  is  used  to  fit  the  density  data  at  Zq  (in  a  least-squares  sense)  with 
the  following  formula: 


n  (z,  L)  =  n  (z.)  exp 


T  1  T  M*o  (r_z0)  (RF  *  z0)  1 

G  (L)  -  1  exp  -  — - - - - - - —  (24) 

L  J  L  R  TJL)  (R-  ♦  z)  J 
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number  density  of  a  particular  gas,  cm" 
altitude,  km 
reference  altitude 

variable  representing  the  geographical  and  geophysical  parameters 
upon  which  density  is  assumed  to  depend 

spherical  harmonic  expansion  described  below 

molecular  weight  of  a  particular  gas,  g/mole 

polar  radius  of  the  earth,  6357  km 

acceleration  of  gravity  at  the  reference  altitude 

-3  2  2 

gas  constant,  8.314  x  10  gm-km  /mole  sec  deg 
exospheric  temperature,  °K. 
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The  function  G(L)  contains  the  information  on  the  variation  of  density  with  selected 
geographic  and  solar  activity  parameters.  It  is  written  as  a  series  of  spherical  harmonics 
whose  coefficients  for  each  gas  have  been  determined  using  the  mass  spectrometer  data 
and  are  tabulated  in  the  model.  The  spherical  harmonics  represent  the  following  kinds  of 
variations  of  density  in  the  thermosphere:  variations  with  solar  activity  (F1Q  7), 
variations  with  geomagnetic  activity  (A^),  variations  with  time  of  the  year  (including 
annual  and  semiannual  terms),  variations  with  local  time  and  latitude  (including  diurnal, 
semidiurnal  and  terdiurnal  terms).  A  total  of  37  coefficients  for  each  thermospheric 
constituent  is  involved.  The  input  parameters  required  by  the  model  are  latitude,  local 
time  and  altitude  of  the  point  in  question,  day  of  the  year,  the  values  of  Fjq  7  both  one 
day  before  the  day  in  question  and  averaged  over  three  solar  rotations,  and  the 
geomagnetic  activity  index  A^,  with  a  six  hour  lag  in  the  OGO-6  model,  and  a  latitude- 
dependent  lag  of  4.8-  :$i/45°  ,  where  0  is  the  geographic  latitude,  in  the  ESRO-4 
model. 

The  merit  of  a  model  based  on  a  spherical  harmonic  analysis  is  that  these 
functions  form  a  complete  set  that  can,  in  principle,  represent  any  degree  of  complexity 
in  the  data  by  a  systematic  increase  in  the  number  of  terms  used.  Furthermore,  since 
the  spherical  harmonics  are  approximate  eigenfunctions  of  the  thermosphere,  relatively 
few  terms  should  be  required.  The  form  of  the  dependence  of  density  on  FjQ  ?  is  similar 
to  that  of  the  Jacchia  1971  model,  except  that  allowance  is  made  for  a  possible  non¬ 
linearity  in  the  correlation  with  the  27-day  F10  ?  variation;  such  an  effect  has  been 
suggested  by  some  data  from  ground-based  observations.  The  dependence  of  density  on 
Ap  included  in  these  models  is  small;  the  possibility  of  a  latitude  variation  of  the 
geomagnetic  activity  effect  has  been  included. 


Once  the  density  distributions  at  the  reference  altitude  and  T®  are  determined 
using  equation  (24),  the  following  formulas  are  used  to  evaluate  the  densities  and  temper¬ 
atures  at  all  other  altitudes: 


T(z,  L)  *  Td  (L)  -  [t^  (L)  -T1m|  exp  [-s  (z-120)] 


(25) 


n(z,  L)  =  n  (120)  D  (z,  T  (L),  s) 

CD 


(26) 


D  (z,  T  (L),  s) 


1  ♦  a  ♦  y 

exp  (-  a  y  £  ) 


(27) 


where  Tj2q  =  temperature  at  120  km  altitude,  355°K 


s  =  temperature  gradient  parameter,  0.022  km 
a  =  s  ♦  1.5  x  10'4 

a  *  l-T120/T«(L) 

K  =  (z  -  120)  (Re  ♦  120)/(Rp  ♦  z) 

y  .  Mg1I0/  «RT.(U 

a  =  thermal  diffusion  coefficient. 
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The  expression  for  T  (z,L)  leads  to  a  temperature  which  increases  with  increasing 
altitude  above  the  120  km  level,  asymptotically  approaching  the  exospheric  temperature 
Ts(L)  at  great  heights.  The  formulations  of  equations  (26)  and  (27)  are  based  on  the 
assumption  of  hydrostatic  equilibrium  above  the  120  km  level,  where  variations  in  all 
quantities  other  than  the  temperature  and  the  Ng  density  are  allowed. 

Accuracy.  The  accuracy  of  these  models  is  on  the  order  of  the  experimental  error 
for  He  and  O  (20%  for  O  and  larger  for  He)  and  about  three  times  the  experimental  error 
for  Nj.  The  root-mean-square  error  of  the  models  (with  respect  to  the  data  which  were 
used  in  their  generation)  is  on  the  order  of  15-20%. 

Limitations.  There  are  two  aspects  of  the  mass  spectrometer  models  which  limit 
their  usefulness.  These  limitations  are  basically  due  to  characteristics  of  the  data  bases 
used  in  their  generation.  Since  each  model  used  data  from  only  about  two  years  of 
satellite  operation,  the  coverage  in  latitude  and  local  time  at  specified  values  of  Pjq  7 
and  Ap  is  not  really  adequate  for  producing  a  definitive  model.  (This  is  especially 
evident  with  respect  to  Ap,  leading  to  a  restriction  of  the  models  to  quiet  geomagnetic 
conditions.)  The  second  limitation  is  related  to  the  coverage  in  altitude.  Since  the  data 
used  in  the  models  were  taken  above  a  particular  altitude,  no  actual  measurements  of 
density  below  these  altitudes  could  be  used.  The  model  predictions  at  lower  altitudes  are 
thus  extrapolations  of  the  model  under  the  assumption  of  hydrostatic  equilibrium. 
Hydrostatic  equilibrium  for  an  individual  gas  species  probably  does  not  occur  below 
~200  km  altitude,  except  in  the  case  of  molecular  nitrogen.  Thus,  these  models  may  be 
suspect  at  all  altitudes  below  those  at  which  data  were  taken,  and  certainly  are 
unreliable  below  -200  km. 
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The  Mass  Spectrometer-Incoherent  Scatter  (MS1S)  Model 


Introduction.  The  OGO-6  and  ESRO-4  models  were  the  first  empirical  models  of 
the  upper  atmosphere  to  be  based  on  data  taken  by  satellite-borne  mass  spectrometers. 
These  models  suffer  from  restricted  data  bases,  however,  and  are  useful  only  when  values 
of  atmospheric  properties  are  required  for  the  ranges  of  altitude,  solar  and  geomagnetic 
activity,  etc.,  which  were  included  in  these  data  bases. 

Since  the  launch  of  these  two  satellites,  a  number  of  additional  spacecraft 
carrying  mass  spectrometers  have  been  placed  into  earth  orbit.  These  satellites,  taken  in 
combination,  provide  a  significant  improvement  in  data  coverage.  For  example, 
measurements  of  atmospheric  density  and  composition  are  now  available  over  a  wide 
range  of  altitudes  (  M40-600  km),  and  for  a  range  of  FjQ  ?  of  ~75-180  x  10  W/m  Hz. 
Furthermore,  a  large  data  base  of  upper-atmospheric  temperature  measurements  has 
been  compiled  at  several  ground-based  incoherent  scatter  radar  sites;  these  data  provide 
complementary  information  in  the  form  of  direct  measurements  of  temperature  with 
good  local  time  coverage  at  particular  latitudes  over  a  wide  range  of  solar  activity  and 
all  seasons. 

The  mass  spectrometer-incoherent  scatter  (MSIS)  model  (Hedin  et  al.,  1977a, b) 
makes  use  of  measurements  of  upper-atmospheric  composition  from  mass  spectrometers 
on  five  satellites  (AE-B,  OGO-6,  San  Marco  3,  Aeros  A  and  AE-C)  and  neutral  tempera¬ 
tures  inferred  from  incoherent  scatter  measurements  at  four  ground  stations  (Arecibo, 
Puerto  Rico;  Jicamarca,  Peru;  Millstone  Hill,  Mass.;  and  St.  Santin,  France).  The  overall 
data  set  covers  the  time  period  from  the  end  of  1965  to  mid-1975  (effectively  a  complete 
11  year  solar  cycle),  values  of  Fj0  7  of  75-180,  values  of  Ap  up  to  -^100,  and  is  moat 
applicable  to  the  altitude  range  200-600  km. 


Model  Formulation.  The  model  formulation  is  virtually  identical  to  that  of  the 
OGO-6  and  ESRO-4  models.  The  following  exceptions  allow  more  versatility  in  the  MSIS 
model: 

1)  The  temperature  and  density  at  120  km  are  functions  of  L,  the  variable  repre¬ 
senting  the  geographical  and  geophysical  input  parameters. 

2)  The  temperature  gradient  parameter  s  is  a  function  of  L. 

3)  The  variation  with  daily  geomagnetic  index  Ap  is  nonlinear,  allowing  the  model  to 
be  applied  to  moderately  disturbed  times. 

4)  A  time-independent  North-South  asymmetry  is  allowed  in  the  densities. 

5)  Because  of  the  likelihood  of  noticeable  departures  from  hydrostatic  equilibrium  in 
the  diurnal  and  magnetic  activity  components  of  the  density  below  200  km,  some 
of  the  terms  in  the  function  G(L)  are  multiplied  by  a  function  of  altitude  which 
approximates  the  effects  of  a  vertical  diffusion  velocity. 

6)  Many  parameters  in  the  function  G(L)  as  formulated  in  the  OGO-6  and  ESRO-4 
models  were  found  to  be  of  little  significance  and  are  set  to  zero  in  the  MSIS 
model. 

Accuracy.  The  uncertainty  in  the  absolute  average  exospheric  temperature  in  the 
model  is  about  15-20°K.  The  coefficient  of  temperature  variation  with  FJ0  7  is  3% 
higher  than  in  the  Jacchia  1971  model  and  10%  higher  than  that  in  the  OGO-6  model. 
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The  MSIS  model  N2  density  at  120  km  is  25%  below  that  in  the  Jacchia  1971  model.  The 
average  temperature  gradient  parameter  s  =  0.0255  gives  a  temperature  gradient  at 
120  km  of  16.7°K/km,  which  is  higher  than  the  10°K/km  of  the  Jacchia  1971  model,  but 
which  agrees  with  the  incoherent  scatter  measurements. 

The  various  atmospheric  species  can  be  crudely  ranked  as  to  the  absolute  error  in 
the  model  representation  of  the  mean  density  profile  above  150  km:  N2(~15%), 
O  ( ~20%),  He  (  — 20%),  Ar(~25%),  H(~50%)  and  O2(~50%).  Fortunately,  the  atmo¬ 
spheric  gas  above  150  km  is  composed  primarily  of  those  species  whose  errors  are  the 
lowest. 

At  120  km,  the  total  density  is  30%  below  that  of  the  Jacchia  1971  model,  pri¬ 
marily  because  of  the  lower  N2  density,  which  is  closer  to  that  measured  by  rocket-borne 
mass  spectrometers  than  is  that  in  the  Jacchia  1971  model.  At  altitudes  above  205  km, 
the  MSIS  model  average  density  is  about  15%  above  that  of  the  Jacchia  1971  model,  well 
within  errors  of  absolute  calibration  and  analysis  for  both  mass  spectrometers  and  orbital 
drag  studies. 

Limitations.  The  MSIS  model  is  believed  to  provide  a  better  representation  of 
global  temperature  and  composition  variations  than  do  other  existing  global  models.  The 
limitations  of  the  model  arise  primarily  from  incomplete  data  at  low  altitudes  and  during 
times  of  disturbed  magnetic  conditions. 
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Jacchia  1977 


The  Jacchia  1977  model  is  the  latest  in  a  series  of  empirical  models  of  the 
thermosphere.  The  previous  (1971)  Jacchia  model  enjoyed  widespread  use,  but  the 
analysis  of  mass  spectrometer  data  from  OGO-6  revealed  that  revision  was  necessary. 
Analysis  of  these  data  showed  that  the  phase  of  the  diurnal  variation  was  a  function  of 
altitude  and,  more  important,  at  a  fixed  altitude  the  phases  of  the  diurnal  variations  of 
the  individual  gases  were  different.  A  major  goal  of  the  model  was  to  reproduce  the 
relative  concentrations  of  N2  and  O  at  450  km  as  derived  from  OGO-6  while  agreeing 
with  the  total  density  determined  by  satellite  drag. 

Jacchia  1977  is  really  two  sets  of  models,  the  static  models  and  a  model  which 
includes  known  thermospheric  variations.  The  static  models  are  indexed  by  exospheric 
temperature,  from  which  a  temperature  profile  is  determined.  The  profile  starts  with  a 
constant  value  of  188°K  at  90  km,  and  an  altitude  gradient  of  zero,  rises  to  an  inflection 
point  at  a  fixed  altitude  of  125  km,  and  then  continues  to  rise,  asymptotically  ap¬ 
proaching  the  selected  exospheric  temperature  at  high  altitudes.  In  addition  to  temper¬ 
ature,  the  static  models  give  number  densities  as  a  function  of  altitude  for  Nj,  Oj,  Ar, 
O,  He,  and  H.  A  fraction  by  volume  for  each  gas  at  90  km  which  agrees  with  the  sea 
level  value  (except  for  a  correction  to  the  helium  fraction  which  is  required  in  order  to 
obtain  proper  thermospheric  helium  results)  is  specified.  Between  90  and  100  km  a 
temporary  empirically  determined  mean  molecular  weight  profile  is  given  and  used,  along 
with  the  temperature  profile,  to  integrate  the  barometic  equation  for  total  mass  density. 
The  departure  from  constant  mean  molecular  weight  is  ascribed  entirely  to  the 


dissociation  of  02;  the  relative  number  densities  of  N^,  Ar,  and  He  are  held  constant  at 
the  90  km  values.  Small  corrections  are  then  applied  to  the  O  and  profiles  in  order 
to  bring  them  into  better  agreement  with  experimental  measurements.  The  final  mean 
molecular  weight  profile  is  computed  from  the  number  densities  thus  determined.  Above 
100  km  the  individual  number  densities  decrease  with  altitude  at  a  rate  determined  by  a 
diffusion  equation  which  includes  thermal  diffusion  (He  and  H)  and  vertical  flux  (H  only) 
terms  as  well  as  the  usual  molecular  diffusion  and  gravity  terms.  For  H  the  boundary 
condition  is  specified  at  500  km  and  is  a  function  of  exospheric  temperature,  unlike  the 
other  gases  for  which  the  boundary  conditions  are  determined  at  100  km  as  described 
above. 

The  Jacchia  1977  model  which  includes  thermospheric  variations  is  an  elaboration 
of  the  static  models  and  attempts  to  model  several  types  of  variation: 

1.  Variation  with  the  solar  cycle 

2.  Variation  with  the  daily  change  in  activity  on  the  visible  disk  of  the  sun 

3.  The  daily,  or  diurnal,  variation 

4.  Variation  with  geomagnetic  activity 

5.  Seasonal-latitudinal  variations 

6.  The  semi-annual  variation. 

The  discussion  of  the  model  in  this  report  is  necessarily  brief,  and  the  complete 
description  (Jacchia,  1977)  should  be  consulted  for  detaiLs. 
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The  mean  exospheric  temperature  is  determined  by  a  function  of  the  daily  F10.7 
solar  flux  (evaluated  at  a  time  delayed  an  amount  depending  on  local  time)  and  a 
Gaussian  weighted  average  of  FjQ  ?  (using  a  time  constant  of  71  days). 

The  diurnal  variation  of  the  temperature  about  its  mean  value  is  given  by  a 
function  of  solar  declination,  latitude,  and  local  time.  For  use  in  determining  number 
density  profiles,  a  set  of  "pseudo-temperature"  profiles  are  introduced,  one  for  each 
constituent.  The  equations  which  describe  these  pseudo-temperature  curves  all  have  the 
same  form,  but  the  phase  of  the  local  time  variation  is  different  for  the  pseudo- 
temperature  associated  with  each  of  the  individual  constituents.  This  is  done  so  that, 
when  integrating  the  diffusion  equation  above  100  km,  the  times  of  peak  density  of  the 
various  constituents  can  be  at  different  hours  of  the  day.  The  phases  of  the  various 
pseudo-temperatures  arc  related  in  such  a  way  as  to  result  in  the  peak  in  total  density 
always  occurring  at  the  same  local  time,  independent  of  altitude. 

Variations  with  geomagnetic  activity  (characterized  by  the  planetary  index  Kp) 
are  divided  into  three  categories.  The  first  category,  the  "thermal  component,"  is 
accounted  for  by  an  increase  in  temperature  at  all  altitudes.  The  amount  of  atmospheric 
heating  which  causes  this  temperature  increase  is  characterized  by  Kp  evaluated  at  a 
time  which  allows  for  a  delay  ckie  to  the  propagation  time  of  the  heating  from  the  polar 
regions  to  the  point  of  interest.  This  delay  time  increases  with  decreasing  invariant 
magnetic  latitude.  A  correction  to  the  temperature  profile  at  all  altitudes  is  determined 
by  the  delayed  Kp  value  and  the  invariant  magnetic  latitude.  This  corrected  temperature 
profile,  of  course,  is  used  when  integrating  the  barometric  and  diffusion  equations.  It 
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might  be  noted  that  this  change  greatly  complicates  use  of  the  model  because  the 
resulting  temperature  profile  is  not  the  same  as  the  static  model  temperature  profile 
with  the  same  exospheric  temperature.  This  means  that  for  Kp  *  0  it  is  not  sufficient  to 
calculate  the  exospheric  temperature,  consult  the  static  model  tables  corresponding  to 
that  exospheric  temperature,  and  then  modify  the  tabulated  densities  as  prescribed  by 
the  model. 

The  second  component  of  the  variations  Aie  to  geomagnetic  activity  is  the  "effect 
of  a  change  in  the  height  of  the  homopause."  This  is  accounted  for  by  a  correction  to  the 
individual  number  densities  which  increases  the  abundance  of  Ar  and  02  while  decreasing 
the  abundance  of  O  and  He;  no  correction  is  made  to  N2  or  H. 

The  third  component  is  the  "equatorial  wave"  which  describes  the  "piling  up"  of 
the  atmosphere  near  the  equator  due  to  meridional  flow  driven  by  the  polar  region 
heating.  ITiis  latitude-dependent  correction  increases  all  constituents  by  the  same 
fractional  amount. 

The  seasonal-latitudinal  variations  are  divided  into  a  "thermospheric"  and  a 
"mesospheric"  variation.  The  thermospheric  term  adjusts  each  gas  by  an  amount 
proportional  to  the  product  of  solar  declination  and  sine  of  the  latitude  of  the  point  of 
interest. 

Although  in  the  troposphere  and  stratosphere,  as  well  as  the  upper  thermosphere, 
the  temperature  is  warmer  in  summer  than  winter,  in  the  mesosphere  the  temperature  is 
warmer  in  winter  than  in  summer.  A  corresponding  behavior  (with  opposite  phase)  is 
expected  for  the  density.  The  mesospheric  seasonal-latitudinal  correction  is  of  no 
importance  above  170  km. 
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The  semi-annual  variation  is  accounted  for  in  the  same  manner  as  in  the  Jacchia 
1971  model  where  the  density  of  all  constituents  is  modified  by  a  percent  which  is  the 
product  of  a  function  of  altitude  and  a  function  of  day  of  the  year. 

Limitations  of  the  Model.  The  great  majority  of  the  data  which  have  been  used  to 
develop  the  model  is  from  altitudes  above  250  km.  The  model  may  therefore  not 
represent  variations  in  the  lower  thermosphere  as  well  as  variations  in  the  higher  altitude 
regions.  The  model  is  not  well  suited  for  use  without  a  computer  because  of  its  com¬ 
plexity.  Relative  to  other  recent  models,  Jacchia  1977  is  not  well  suited  for  use  on  a 
computer  either.  The  modification  of  the  basic  form  of  the  temperature  profile  because 
of  geomagnetic  activity  requires  either  numerical  integration  of  the  barometric  and 
diffusion  equations  for  each  data  point  of  interest,  or  storage  of  a  vast  amount  of  pre¬ 
computed  data.  The  introduction  of  psuedo-temporatures  and  the  later  correction  of 
densities  destroys  much  of  the  aesthetic  value  which  might  have  been  derived  from 
basing  density  profiles  on  the  diffusion  equation.  One  might  well  ask  whether  the  great 
complexity  relative  to  other  modeLs  is  adequately  justified  by  any  compensating  sig¬ 
nificant  incease  in  accuracy.  Computational  ease  and  efficiency  have  apparently  not 
been  given  much  weight  in  the  development  of  the  model. 

Accuracy.  No  adequate  discussion  of  the  accuracy  of  the  model  is  presented  by 
its  author,  but  comparison  with  experimental  data  (see  the  next  section)  leads  to  an 
estimated  r.m.s.  error  for  densities  of  about  19%. 
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IIL  DISCUSSION  AND  COMPARISON  OF  MODELS 


Two  major  factors  which  determine  the  usefulness  of  a  model  are  its  convenience 
and  its  accuracy.  Although  the  need  for  accuracy  is  obvious,  the  need  for  convenience  is 
often  not  clearly  recognized.  To  a  large  extent,  the  requirements  for  a  model  to  be 
convenient  for  use  are  the  same  whether  calculations  are  to  be  made  manually  or  by  an 
electronic  computer.  A  notable  exception  is  in  the  use  of  extensive  tables,  which  are 
more  appropriate  for  the  human  than  the  computer.  A  classic  example  is  the  evaluation 
of  an  elementary  trigonometric  function  which  a  human  might  look  up  in  a  table  of  pre¬ 
calculated  values;  the  storage  space  which  would  be  taken  up  by  such  a  table  in  a 
computer  is  too  valuable  for  such  use,  and  the  elementary  functions  are  therefore 
calculated  from  a  formula.  Hereafter  we  will  assess  convenience  in  terms  of  imple¬ 
mentation  of  the  model  for  use  by  a  computer. 

A  model  should  have  a  simple  structure.  It  can  then  be  implemented  with  a 
compact  code  which  requires  little  storage  space.  This  desire  for  simple  structure  is 
somewhat  in  conflict  with  the  complicated  behavior  of  the  thermosphere  which  tends  to 
require  more  complexity  in  the  model  if  adequate  accuracy  is  to  be  maintained.  Simple 
structure  also  usual'y  leads  to  fast  execution  time;  however  a  structure  which  is 
conceptually  simple  can  lead  to  slow  execution  in  some  cases.  The  most  common 
instance  of  this  difficulty  is  in  the  semi-empirical  models  which  specify  a  temperature 
profile  and  then  give  a  differential  equation  (involving  temperature)  which  must  be 
solved  to  obtain  density  profiles.  The  differential  equation  is  solved  by  numerical 
integration,  a  simple  but  time-consuming  process.  If  the  equation  can  be  integrated 
analytically  (as  is  the  case  when  using  a  "Bates-Walker"  approach),  the  usefulness  of  the 
model  is  improved. 
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The  model  must  provide  the  required  parameters.  Moat  modern  models  provide 
temperature  and  densities  of  the  major  species.  If  only  total  density  is  required,  this 
detailed  output  may  have  been  purchased  at  some  cost  in  terms  of  efficiency.  The  input 
parameters  required  by  the  model  must  be  available.  Such  indices  as  Fjg  ^  and  have 
been  used  in  spite  of  their  well-known  inadequacies  largely  because  of  their  availability. 

That  a  model  must  be  accurate  is  obvious.  Less  obvious  is  how  the  accuracy  is  to 
be  characterized  for  a  particular  use.  A  perfect  description  of  the  atmosphere  is  the 
ideal,  but  some  kinds  of  inaccuracies  are  unimportant  for  specific  applications.  A 
satellite  system  designer  who  is  required  to  estimate  the  orbital  lifetime  until  reentry  of 
a  satellite  designed  for  a  two  year  lifetime  may  be  quite  satisfied  with  a  model  having 
the  correct  mean  value  when  averaged  over  a  month's  time  span  even  if  the  errors 
averaged  over  a  single  orbit  are  quite  large.  On  the  other  hand,  a  tracking  station 
operator  may  attach  no  importance  at  all  to  the  long-term  mean  behavior  of  the 
atmosphere  (which  can  be  empirically  adjusted  for  in  day-to-day  operations)  if  the  model 
could  only  accurately  predict  the  short-term  variations. 

Now  we  will  try  to  make  some  assessment  concerning  the  usefulness  of  the  models 
presented  in  this  report.  Models  pnxfciced  prior  to  1970  were  developed  from  data  bases 
with  significant  limitations.  Most  of  the  data  were  from  satellites  with  perigee  altitudes 
above  250  km  and  were  based  on  drag  calculations  which  only  give  total  density.  The 
only  advantage  these  models  have  over  more  recent  models  is  their  greater  simplicity 
and  the  many  years  of  experience  which  some  groups  have  had  in  using  them.  More 
recent  models  are  more  suitable  except  for  those  situations  where  only  total  density  is 
required  and  some  loss  of  accuracy  can  be  tolerated  in  exchange  for  small,  fast  computer 
code  or  where  the  cost  of  replacing  present  software  greatly  exceeds  the  benefits 
achievable  by  the  use  of  better  (but  still  far  from  perfect)  models. 


The  OGO-6  and  ESRO-4  models  have  the  same  structure  as  the  MSIS  model  but 
were  based  on  more  limited  data,  and  these  two  models  can  probably  be  regarded  as 
obsolete. 

The  McDonnell  Douglas  model  (MDAC)  has  as  its  most  interesting  and  potentially 
important  feature  a  specific  incorporation  of  high  latitude  heating  effects.  The  present 
version  of  the  model,  unfortunately,  does  not  specify  how  the  parameters  which  describe 
this  heating  are  to  be  derived  from  available  geophysical  data,  and  only  "typical"  values 
for  these  parameters  are  given.  Thus,  the  model's  greatest  potential  asset  is  presently 
one  of  its  greatest  weaknesses,  and  until  the  model  is  developed  further  it  cannot  be 
recommended  for  general  use. 

The  best  models  for  general  use  are  probably  Jacchia  1971  (or  the  nearly  identical 
CIRA  1972)  and  MSIS.  Jacchia  1971/CIRA  1972  is  an  internationally  recognized  model 
which  is  familiar  to  most  workers  in  the  field.  It  is  thoroughly  documented,  and  an 
excellent  FORTRAN  version  of  the  computer  code  has  been  published  in  the  CIRA  1972 
report.  It  is  primarily  based  on  <frag  data  from  satellites  with  perigee  altitudes  above 
200  km.  It  does  not  incorporate  some  of  the  most  recent  knowledge  of  the  behavior  of 
individual  species  and  thus  is  especially  of  value  when  only  total  density  at  altitudes 
above  200  km  is  required.  The  model  requires  numerical  integration  of  the  diffusion 
equation  which  is  a  time-consuming  procedure;  in  order  to  avoid  this  difficulty  some 
users  have  developed  modified  versions  of  the  model  which  give  approximately  the  same 
results  but  are  computationally  more  efficient. 

Jacchia's  1977  model,  although  still  limited  in  the  amount  of  low  altitude  data  on 
which  it  is  based,  does  describe  the  diurnal  behavior  of  the  individual  species  in  a  manner 


consistent  with  recent  observations.  The  major  objection  to  this  model  is  its  complexity 
and  resulting  inconvenience  of  use.  It  appears  that  development  of  the  model  was  guided 
almost  entirely  by  considerations  of  accuracy. 

The  data  base  used  to  develop  the  MSIS  model  includes  low  altitude  data  and  mass 
spectrometer  determined  measurements  of  individual  species.  Numerical  integration  is 
not  required,  resulting  in  a  reasonably  fast  computer  code.  There  is  no  published  version 
of  the  computer  code.  The  Goddard  group  has  provided  a  FORTRAN  version  of  the 
model  upon  request,  but  there  are  some  difficulties  with  this  code.  It  appears  to  have 
been  developed  over  a  period  of  time  and  is  not  as  "clean"  or  well  written  as  the  CIRA 
1972  code.  It  is  written  in  a  non-ANSI  dialect  of  FORTRAN  which  is  not  acceptable  at 
some  computer  installations  without  non-trivial  modifications.  Despite  these  limitations, 
it  is  probably  the  best  presently  available  model  for  general  use,  especially  for  use  at 
altitudes  below  200  km  and  where  individual  species  data  are  required.  As  the  model  is 
based  on  only  a  fraction  of  the  Atmosphere  Explorer  data,  publication  of  a  revised 
version  in  the  near  future  would  not  be  surprising.  If  we  may  be  guided  by  memory  of  the 
replacement  of  the  OGO-6  model  by  the  MSIS  model,  we  may  anticipate  that  a  revised 
version  of  the  MSIS  model  would  have  a  structure  nearly  identical  to  the  present  version. 

A  specific  comparison  of  a  number  of  models  with  measured  data  is  given  in 
Tables  1-13.  The  data  consist  of  in-situ  measurements  of  total  density  made  by  cold 
cathode  ion  gauges  carried  on  board  NASA's  AE-C,  D,  and  E  satellites.  Data  were  from 
December  1973  to  November  1974  (AE-C),  October  1975  to  January  1976  (AE-D),  and 
November  1975  to  October  1976  (AE-E).  Measurements  were  made  in  an  altitude  range 
of  approximately  135  to  300  km.  For  every  experimental  determination  of  density,  each 
of  the  models  was  used  to  calculate  a  density  using  values  of  geophysical  and  solar 
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parameters  appropriate  for  the  time  and  location  of  the  experimental  measurement.  The 
quantity 

r  =  log  /  measured  density  \ 
e  \  model  density  / 

was  then  used  to  compute  three  statistics  which  may  be  used  to  examine  the  ability  of 
the  models  to  predict  thermospheric  density.  Let  us  denote  by  r^  this  quantity  calculated 
for  the  i-th  density  determination  in  the  data  base. 

The  mean,  or  average,  error  of  a  model  is  characterized  by 

n 

T  1  H  E  ri  • 

1=1 


where  n  is  the  number  of  {joints  used.  The  r.m.s.  error  is  characterized  by 


and  the  r.m.s.  error  after  removing  the  mean  error  is  characterized  by 


n 


L 

i=i 
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TABLE  1 


Comparison  of  Measured  Densities  to  Model  Values 


All  altitudes  (135-300  km) 
n  =  32,  720 


.  / 

r 

vr 

’DENSEL 

0.152 

0.283 

J65 

-0.186 

0.258 

J71 

0.058 

0.173 

MDAC 

0.230 

0.307 

MSIS 

0.001 

0.163 

J77 

0.054 

0.182 

r  =  loKe 

/ measured  density  \ 
l  model  density  / 

TABLE  2 


Comparison  of  Measured  Densities  to  Model  Values 


135  km 

n  =  331 

/  1 

l~i  -2 

f 

vr 

Vr  “r 

'DENSEL 

0.236 

0.289 

0.166 

J65 

-0.089 

0.191 

0.169 

J71 

-0.040 

0.150 

0.145 

MDAC 

-0.017 

0.150 

0.149 

MS1S 

0.110 

0.170 

0.130 

J77 

-0.052 

0.199 

0.192 

_  _ _ /  measured  density  \ 

'  l0ge  (  modelde"nsTty  ) 


TABLE  4 


Comparison  of  Measured  Densities  to  Model  Value 


150  km 

n  =  1655 

/“* 

i~1  -2 

r 

vr 

Vr-r 

’DENSEL 

0.171 

0.234 

0.159 

J65 

-0.108 

0.191 

0.158 

J71 

0.034 

0.143 

0.139 

MDAC 

0.108 

0.185 

0.150 

MSIS 

0.057 

0.135 

0.122 

J77 

0.031 

0.165 

0.162 

r 


l0*«( 


measured  density  \ 
model  density  } 
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TABLE  5 


Comparison  of  Measured  Densities  to  Model  Value 


160  km 
n  =  2059 


Fi 

r 

Vr 

% 

V r  •+ 

'DENSEL 

0.178 

0.238 

0.159 

J65 

-0.132 

0.200 

0.150 

J71 

0.057 

0.152 

0.141 

MDAC 

0.164 

0.231 

0.163 

MSIS 

0.041 

0.130 

0.123 

J77 

0.050 

0.168 

0.161 

,  _  /measured  density  1 
,ogc  (  model  d^sitT  I 


r 


TABLE  6 


Comparison  of  Measured  Density  to  Model  Values 


170  km 

n  =  1371 

r T 

jl 

ii 

KJ 

r 

Vr 

Vr" 

'DENSEL 

0.223 

0.275 

0.161 

J6S 

-0.133 

0.200 

0.149 

J71 

0.082 

0.169 

0.148 

MDAC 

0.216 

0.276 

0.172 

MSIS 

0.037 

0.140 

0.135 

J77 

0.069 

0.166 

0.151 

/  measured  density  \ 
"  lo*e  (  mS^fr'nsity  ) 
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TABLE  7 


Comparison  of  Measured  Densities  to  Model  Values 


180  km 

n  =  1367 

[~i 

,  /  5T* 

r 

vr 

V 

•DENSEL 

0.256 

0.308 

0.172 

J6S 

-0.149 

0.214 

0.154 

J71 

0.091 

0.178 

0.152 

MDAC 

0.253 

0.310 

0.178 

MS1S 

0.029 

0.145 

0.142 

J77 

0.076 

0.175 

0.158 

r 


/  measured  density  \ 
'°*e  l  mo<diei  densfty  • 
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TABLE  8 


Comparison  of  Measured  Densities  to  Model  Values 


200  km 

n  =  1982 

r 

/ 

r 

>/** 

’DENSEL 

0.286 

0.340 

0.184 

J65 

-0.188 

0.245 

0.156 

J71 

0.099 

0.185 

0.156 

MDAC 

0.312 

0.365 

0.188 

MSIS 

0.025 

0.152 

0.150 

J77 

0.091 

0.189 

0.165 

*  1oM 


,  measured  deasity  | 
model  density 


; 
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TABLE  9 


Comparison  of  Measured  Densities  to  Model  Values 


220  km 

n  =  1100 

/  ^ 

r 

V  r 

'DF.NSEL 

0.084 

0.238 

J6S 

-0.266 

0.324 

J71 

0.039 

0.185 

MDAC 

0.277 

0.350 

MSIS 

-0.068 

0.185 

J77 

0.033 

0.188 

r 


,  measured  density  i 
(  model  density  ' 


TABLE  10 


Comparison  of  Measured  Densities  to  Model  Values 


240  km 

n  =  804 

/jt 

.  /  2  2 

r 

vr 

'DENSEL 

-0.108 

0.251 

0.226 

J65 

-0.332 

0.375 

0.174 

J71 

-0.006 

0.182 

0.182 

MDAC 

0.2S2 

0.332 

0.215 

MSIS 

-0.145 

0.226 

0.174 

J77 

-0.009 

0.178 

0.178 

r 


, _  /  measured  density  \ 

IO*e  (  ■'  model density  > 
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TABLE  11 

Comparison  of  Measured  Densities  to  Model  Value 


260  km 

n  =  675 

/  2 

,  17  -2 

r 

Vr 

V,-r 

'DFNSEL 

-0.213 

0.300 

0.212 

J65 

-0.316 

0.358 

0.169 

J71 

0.023 

0.182 

0.181 

MDAC 

0.290 

0.360 

0.214 

MSIS 

-0.141 

0.224 

0.174 

J77 

0.035 

0.184 

0.180 
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TABLE  12 


Comparison  of  Measured  Densities  to  Model  Values 


280  km 

n  =  625 

r 

'DENSEL 

-0.262 

0.337 

0.211 

J65 

-0.283 

0.331 

0.172 

J71 

0.068 

0.192 

0.179 

MDAC 

0.337 

0.399 

0.215 

MSIS 

-0.112 

0.206 

0.172 

J77 

0.091 

0.200 

0.178 

r 


•  I  measured  density  > 
^e'  model  density  ' 
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Perhaps  the  most  notable  (and  most  discouraging)  feature  of  Table  1  is  the  lack  of 
substantial  improvement  in  model  development  since  1971.  In  fact,  except  for  an  overall 
multiplicative  factor  due  to  being  based  on  data  taken  at  times  near  maximum  solar 
activity,  one  might  well  conclude  that  there  has  been  little  improvement  since  1965.  If 
individual  species  were  examined  instead  of  total  density,  the  progress  of  the  last  several 
years  would  presumably  be  more  apparent.  J71,  MSIS  and  J77  all  give  similar  accuracy 
when  used  to  predict  total  density  and  are  significantly  superior  to  'DENSEL,  J65  and 
MDAC.  For  predicting  number  densities  of  individual  species  it  may  be  anticipated  that 
the  performance  of  J71  would  be  poorer  than  that  of  MSIS  and  J77. 

After  nearly  twenty  years  of  thermospheric  modeling  work,  it  is  interesting  to  ask 
about  the  accuracy  of  the  models.  This  is  a  question  not  often  discussed  in  any 
significant  detail  by  the  authors  of  the  various  modeLs.  An  indication  of  the  limits  of 
recent  modeLs  is  seen  in  Figure  10,  taken  from  the  OGO-6  model  report.  The  comparison 
seen  is  between  the  model  prediction  and  the  data  on  which  the  model  is  based.  If 
compared  with  other  data  it  would  undoubtedly  be  found  that  the  agreement  would  be 
even  worse.  Some  of  the  scatter  in  Figure  10  is  attributable  to  measurement  error,  but 
it  seems  clear  that  much  of  it  is  due  to  limitations  of  the  model.  Figures  such  as  this 
clearly  demonstrate  the  limitations  of  present-day  models.  Although  this  performance 
may  seem  rather  poor  in  view  of  the  ever  greater  accuracy  demanded  by  users,  a  glance 
at  Figure  11  may  serve  as  a  reminder  that  the  remaining  errors  are  not  especially  large 
when  compared  with  the  great  variations  present  in  the  upper  atmosphere.  If  one  could 
predict  the  stock  market  behavior  as  accurately  as  thermospheric  density,  his  future 
financial  security  would  be  well  assured. 
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Figure  10, 


Ratio  of  the  measured  density  to  the  model  Nj  density  as  a  function  of 
time  from  perigee,  geographic  latitude,  local  time,  FJ0  ?  solar  flux, 
magnetic  index  A  and  day  of  year  (from  Hedin  et  al.,  1972). 
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(O-CI 


Figure  It.  Observed  and  computed  densities,  reduced  to  a  standard  height  of  530  km 
for  satellite  1964  76 A  (Explorer  24).  Computed  densities  were  derived 
using  the  Jacchia  1971  model  (from  Jacchia,  1971). 
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Why  is  it  that  models  do  not  do  a  better  job  of  predicting  thermospheric  density? 
Part  of  the  answer  lies  in  the  selection  of  parameters  to  characterize  conditions  in  the 
thermosphere.  For  example,  the  amount  of  EUV  heating  is  often  measured  by  the  10.7- 
cm  solar  flux.  The  10.7-cm  flux  itself  cannot  significantly  heat  the  upper  atmosphere.  It 
has  been  established,  however,  that  there  is  a  reasonable  correlation  between  the  EUV 
flux  and  the  10.7-cm  flux.  This  correlation  is  only  approximate  and  may  not  be  adequate 
for  use  in  accurate  models.  However,  the  10.7-em  flux  is  measured  regularly  and  hence 
is  readily  available;  better  a  "fair"  index  than  none  at  all.  A  similar  situation  exists  with 
respect  to  the  energetic  particle  heating  which  occurs  during  geomagnetic  storms.  This 
heating  is  related  to  disturbances  in  the  earth's  magnetic  field,  but  the  relationship  is 
undoubtedly  not  adequate  for  high  accuracy  models.  There  is  a  need  for  new  parameters 
to  be  made  available  which  more  directly  characterize  phenomena  which  affect  the 
thermosphere.  Although  this  need  has  been  recognized  for  over  a  decade,  no  significant 
new  parameters  have  emerged. 

Even  if  such  improved  parameters  were  available,  it  is  not  clear  that  extremely 
accurate  models  are  possible.  Gravity  waves  and  other  transient  disturbances  affect  the 
density  in  a  seemingly  random  fashion.  It  may  be  that  ultimately  we  will  be  able 
accurately  to  predict  the  average  state  of  the  atmosphere  but  not  be  able  to  predict 
short-term  (e.g.  on  a  time  scale  shorter  than  a  few  hours)  departures  from  this  average 
state.  In  the  language  of  the  meteorologist  it  may  be  that  we  can  predict  the  "climate" 
but  not  the  "weather."  It  is  not  clear  at  present  what  the  ultimate  accuracy  of  models 
can  be.  There  is  still  much  to  be  done  in  developing  the  theoretical  understanding  of 
thermospheric  behavior,  particularly  the  interaction  of  the  thermosphere  with  the  next 
lower  altitude  region,  the  mesosphere. 
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Although  there  is  no  immediate  prospect  of  a  revolutionary  improvement  in  the 
accuracy  of  thermospheric  models,  there  are  nevertheless  good  reasons  for  believing  that 
prediction  of  thermospheric  behavior  may  undergo  significant  improvement  in  the  long 
run.  Unlike  tropospheric  weather  prediction,  there  does  not  exist  the  problem  of  local 
geography  causing  small-scale  local  differences.  Propagation  times  for  disturbances  are 
of  the  order  of  hours  which  helps  to  smooth  out  perturbations  quickly.  The  coupling 
between  the  mesosphere  and  the  thermosphere  and  the  unknown  behavior  of  the 
atmosphere  in  the  region  near  90  km  is  a  problem,  as  is  the  lack  of  routine  monitoring  of 
the  state  of  the  entire  thermosphere  by  the  equivalent  of  present  meteorological 
satellites  which  monitor  the  troposphere  and  stratosphere.  As  remote  sensing  techniques 
are  developed  and  applied  to  the  mesosphere  and  lower  thermosphere,  we  may  reasonably 
expect  that  there  will  be  steady,  though  probably  not  spectacular,  progress  toward  higher 
accuracy  models. 
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IV.  DENSITY  CONTOUR  PLOTS 


Contour  plots  with  isopleths  of  constant  total  mass  density  are  shown  in  this 
section  for  the  most  important  modeLs  in  current  use:  'DENSEL,  Jacchia  1971,  MDAC, 
MSIS,  and  Jacchia  1977.  The  000-6  and  ESRO-4  modeLs  are  assumed  to  have  been 
replaced  by  MSIS.  As  the  emphasis  for  this  report  is  on  the  lower  thermospheric  regions 
where  satellite  (frag  is  important,  the  modeLs  have  been  evaluated  at  140  and  180  km  to 
show  the  low  altitude  behavior,  and  at  400  km  to  give  an  indication  of  the  behavior  at 
higher  altitudes. 

The  contour  lines  are  labeled  with  a  two  digit  code,  and  the  densities 
corresponding  to  the  various  code  numbers  are  given  in  a  table  to  the  right  of  each  plot. 
It  should  be  noted  that  each  plot  is  scaled  independently  and  that  in  two  different  plots, 
contour  lines  with  the  same  label  may  correspond  to  different  densities.  The  primary 
coordinates  of  local  time  (in  hours)  and  latitude  (in  degrees)  are  given  on  the  axes  of  the 
plot  for  the  northern  hemisphere.  The  values  of  various  other  parameters  are  given  in 
the  heading  above  the  plot.  Equinox,  Summer,  and  Winter  denote  the  dates  March  21, 
July  21,  and  December  21,  respectively. 

Although  to  the  greatest  extent  possible,  the  plots  have  been  made  for  identical 
conditions  for  each  of  the  models,  some  accommodation  has  been  made  to  the  varying 
requirements  of  the  modeLs.  The  value  labeled  P  10.7  has  been  used  for  both  FJ0  7  (the 
long  term  average)  and  P]Q  „  (the  daily  value).  There  is  no  density  variation  for 
T)ENSEL  at  140  km,  so  contour  plots  at  that  altitude  are  not  shown.  MSIS  requires 
specification  of  Ap  rather  than  Kp.  A  value  of  Ap  =  7  has  been  used  for  the  plots  labeled 

K  =  2. 

P 
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In  some  eases  both  geomagnetic  and  geographic  latitude  are  required.  The  plots 
given  were  produced  by  setting  both  the  geomagnetic  and  geographic  latitude  to  the 
same  value  corresponding  to  the  indicated  latitude. 

It  is  interesting  to  note  that  the  computer  version  of  the  model  obtained  from 
JacchiaS  group  does  not  apply  the  method  described  in  Jacchia  1977  for  correcting  the 
temperature  profile  for  geomagnetic  activity,  but  rather  uses  the  static  model 
temperature  profile  which  has  the  correct  exospheric  temperature.  At  the  time  of  this 
writing,  JacchiaS  group  did  not  have  a  computer  routine  which  corresponded  exactly  to 
the  model  description  (private  communication).  The  lower  altitudes  are  most  likely  to  be 
affected  by  this  change,  although  it  is  anticipated  that  differences  should  not  be  large. 
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LABORATORY  OPERATIONS 

The  Uboniorr  Operations  at  The  Aerospace  Corporation  I*  conducting 
esperimanlal  nnd  theoretical  Investigations  wtulitr  lof  tKu  evaluation  nnd 
application  o(  aciantiBc  advances  to  nan  military  cone  apis  and  ayatama.  Var- 
aatllity  and  fleutbility  have  been  developed  to  a  high  dag raa  by  the  laboratory 
personnel  in  dealing  with  the  many  probtema  encountered  in  the  nation's  rapidly 
developing  apace  and  miaeile  systems.  Eaperttae  In  the  tateel  scientific  devel¬ 
opments  ta  vital  to  the  accomplishment  a i  leeks  related  to  theee  problems.  The 
laboratories  that  contribute  to  this  research  arai 

Aerophysics  Laboratory;  Launch  and  reentry  aerodynamics,  heat  trans¬ 
fer,  reentry  physics,  diemiial  kinetic s.  structural  mechanics,  night  dynamics, 
atmospheric  pollution,  and  high -power  gas  lasers. 

Chemistry  and  Physics  Laboratory:  Atmospheric  reactions  and  atmos- 
pheric  optics,  oiemical  reactions  in  polluted  atmospheres,  chemical  reactions 
of  escited  species  In  rocket  plumes,  chemical  thermodynamics,  plasma  and 
laser  induced  reactions,  laser  < hemiatry.  propulsion  chemistry,  space  vacuum 
and  radiation  effects  on  materials,  lubrication  and  surface  phenomena,  photo¬ 
sensitive  materials  and  sensors,  high  precisian  laser  ranging,  and  the  appli¬ 
cation  of  physics  and  chemistry  to  problems  of  law  enforcement  and  biomedicine. 

Electronics  Research  Laboratory:  Electromagnetic  theory,  devices,  and 
propagation  phenomena,  including  plasma  electromagnetics;  quantum  electronics 
lasers.  and  electre-epttcs.  rommunicellon  sciences,  applied  electronic  s.  semi¬ 
conducting.  superconducting,  and  crystal  device  physics,  optical  and  acoustical 
imaging,  atmospheric  pollution,  millimeter  wave  and  far-infrared  technology. 

Materials  Sciences  laboratory--  Development  of  new  materials,  metal 
mains  composites  and  new  forms  of  carbon,  tsst  and  evaluation  of  graphite 
and  ceramics  in  reentry,  spacecraft  materials  and  electronic  components  In 
nuclear  weapons  environment,  application  of  fracture  mechanics  to  stress  cor¬ 
rosion  and  fatigue-induced  fractures  in  structural  metals. 

Space  Sciences  Laboratory  Atmospheric  and  Ionospheric  physics,  radla- 
lion  from  tile  atmosphere,  density  and  composition  of  the  atmosphere,  aurorae 
and  airglow.  magnetoepheric  physics,  cosmic  rsys,  generation  and  propagation 
of  plasma  waves  In  the  magnetosphere,  solar  physics,  studies  of  solar  magnetic 
fields,  space  astronomy,  s-ray  astronomy,  the  effects  of  nuclear  euptoalooe. 
magnetic  storms,  and  solar  activity  on  the  earth's  atmosphere,  ionosphere,  and 
magnetosphere,  the  effects  of  optical,  electromagnetic,  and  particulate  radia¬ 
tions  in  space  on  space  systems. 
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